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ABSTRACT 

This is the second in a series of papers in which we derive simultaneous constraints 
on cosmology and X-ray scaling relations using observations of massive, X-ray flux- 
selectcd galaxy clusters. The data set consists of 238 clusters with 0. 1-2.4 keV lu- 
minosities > 2.5 X 10^'*/i7Q^ ergs^^, and incorporates follow-up observations of 94 of 
those clusters using the Chandra X-ray Observatory or ROSAT [11 were observed with 
both). The clusters are drawn from three samples based on the i^OS'^T All-Sky Sur- 
vey: the iZ05'^T' Brightest Cluster Sample (78/37 clusters detected/followed- up), the 
ROSAT-YjSO Flux-Limited X-ray sample (126/25), and the bright sub-sample of the 
Massive Cluster Survey (34/32). Our analysis accounts self-consistently for all selec- 
tion effects, covariances and systematic uncertainties. Here we describe the reduction 
of the follow-up X-ray observations, present results on the cluster scaling relations, and 
discuss their implications. Our constraints on the luminosity-mass and temperature- 
mass relations, measured within rsoo, lead to three important results. First, the data 
support the conclusion that excess heating of the intracluster medium (or a combina- 
tion of heating and condensation of the coldest gas) has altered its thermodynamic 
state from that expected in a simple, gravitationally dominated system; however, this 
excess heat is primarily limited to the central regions of clusters {r < 0.15r5oo). Sec- 
ond, the intrinsic scatter in the center-excised luminosity-mass relation is remarkably 
small, being bounded at the < 10 per cent level in current data; for the hot, massive 
clusters under investigation, this scatter is smaller than in either the temperature- 
mass or Yx-mass relations (10-15 per cent). Third, the evolution with rcdshift of 
the scaling relations is consistent with the predictions of simple, self-similar models 
of gravitational collapse, indicating that the mechanism responsible for heating the 
central regions of clusters was in operation before redshift 0.5 (the limit of our data) 
and that its effects on global cluster properties have not evolved strongly since then. 
Our results provide a new benchmark for comparison with numerical simulations of 
cluster formation and evolution. 

Key words: large-scale structure of Universe - X-rays: galaxies: clusters. 



1 INTRODUCTION 

Establishing the relationship between total mass and ob- 
servable quantities is a crucial step in deriving cosmological 
constraints from the growth of cosmic structure using galaxy 
clusters. Not only can these scaling relations provide useful 
proxies for mass, but they are also fundamentally impor- 
tant in accounting for selection effects such as Eddington 
bias and Malmquist bias. 

The construction of X-ray flux-selected cluster samples 
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out to redshift z = 0.5 and beyond has now enabled investi- 
gations of dark energy using these data (Mantz et al. 2008, 
hereafter M08; Vikhlinin et al. 2009b). However, the task of 
calibrating X-ray scaling relations has become correspond- 
ingly more complex; the evolution with cosmic time of the 
scaling relations and their scatter must be well understood, 
since such evolution can be degenerate with the effects of 
dark energy. The effect on cosmological constraints of sys- 
tematic uncertainties in the scaling relations has been dis- 
cussed in the context of future surveys by, e.g., Sahlen et al. 
(2009). 

Fortunately, as we describe below, it is possible to simul- 
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taneously constrain both the evolution of the scahng rela- 
tions and cosmological parameters, using a flux-limited sam- 
ple of which some clusters have been targeted by detailed, 
follow-up X-ray observations. To distinguish X-ray methods 
from other measures of the growth of cosmic structure, in- 
cluding those using optically selected clusters, we refer to 
the resulting data set as the cluster X-ray luminosity func- 
tion (XLF), although, strictly speaking, it contains a great 
deal more information than the luminosity function alone. 

This is the second of a series of papers in which we 
address these issues. In a companion paper (Mantz et al. 
2009, hereafter Paper I) we describe the statistical meth- 
ods required to simultaneously constrain the scaling rela- 
tions and cosmology, self-consistently accounting for all se- 
lection effects^ and systematic uncertainties, and present the 
cosmological results from our analysis. This paper focuses 
on the reduction of the follow-up X-ray observations, and 
the constraints on the scaling relations from the simulta- 
neous analysis. In Papers III (Rapetti et al. 2009) and IV 
(Mantz, Allen, & Rapetti 2009), we respectively apply our 
analysis to investigations of modified gravity and neutrino 
properties. 

In addition to their utility for cosmological investiga- 
tions, cluster scaling relations are of significant astrophysical 
importance. Of primary interest is the heating mechanism 
that prevents cooling gas in dense cluster cores from con- 
densing into stars and molecular gas at much higher rates 
than are observed (for reviews, see Peterson & Fabian 2006; 
McNamara & Nulsen 2007). The shape and evolution of the 
scaling relations, and specifically any departures from the 
simplest predictions for gravitationally dominated systems, 
can provide information on the physical mechanisms respon- 
sible for averting strong cooling and star formation. 

The details of the cluster sample selection, the follow-up 
observations, and their reduction are discussed in Section 2. 
Sections 3 describes the scaling relation model and Section 4 
our statistical method, which is more comprehensively de- 
tailed in Paper I. In Section 5, we present constraints on the 
scaling relations, and investigate various extensions to the 
simplest model, including possible evolution with redshift 
and asymmetric scatter. Section 6 contains a discussion of 
the infiuence of cool. X-ray bright gas in cluster centers on 
the scaling relations. Implications of our results for the as- 
trophysics of the intracluster medium (ICM) are discussed 
in Section 7. 

Unless otherwise noted, masses, luminosities and dis- 
tances in tables and figures are calculated with respect to a 
reference cosmology, defined to be spatially flat, with Hub- 
ble constant h = //o/100kms~^ Mpc~^ = 0.7, present mean 
matter density with respect to the critical density ilm = 0.3, 
and dark energy in the form of a cosmological constant. We 
adopt the conventional deflnition of cluster radius in terms 
of the critical density of the Universe; thus, ta is the radius 
within which the mean density of the cluster is A times the 
critical density at the redshift of the cluster. 
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Figure 1. Luminosity-redshift distribution of clusters in the BCS 
(blue, filled circles), REFLEX (black, open circles) and Bright 
MACS (red triangles) samples which are above the respective 
flux limits (see text). The adopted minimum luminosity of 2.5 X 
ergs ^ is indicated by the dashed line. Luminosities in 
this plot are calculated from survey fluxes, assuming our reference 
cosmology, and refer to the 0.1-2.4 keV band in the cluster rest 
frame. Error bars are not shown. 



2 DATA 

2.1 Cluster samples 

Our data are drawn from three wide-area cluster sam- 
ples derived from the ROSAT All-Sky Survey (RASS; 
Triimper 1993): the 7? 05^ T Brightest Cluster Sample (BCS; 
Ebeling et al. 1998), the ROSAT-¥jSO Flux-Limited X-ray 
sample (REFLEX; Bohringer et al. 2004), and the bright 
sub-sample of the Massive Cluster Survey (Bright MACS; 
Ebeling, Edge, & Henry 2001; Ebeling et al. 2010). Each 
sample covers a distinct volume of the Universe: the BCS 
covers the northern sky at 2 < 0.3, REFLEX covers the 
southern sky^ at 2: < 0.3, and Bright MACS covers higher 
redshifts, 0.3 < 2 < 0.5, at declinations > —40°. The dis- 
tributions of cluster detections from each sample in redshift 
and luminosity are shown in Figure 1 (here luminosities in 
the rest frame 0. 1-2.4 keV band are calculated straightfor- 
wardly from the survey fluxes, assuming our reference cos- 
mology). 

The flux limits, Fum, defining the samples are respec- 
tively 4.4 X 10"^^ 3.0 X 10"^^ and 2.0 x 10"^^ ergs"^ cm"^ 
in the 0.1-2.4 keV T^OS'^T energy band. In order to restrict 
the data set to massive clusters, for which the scaling rela- 
tions can be described simply, we select only clusters whose 
flux in this band satisfies 



F ^ max [Fmiu{z), F^^] , 
where 



(1) 



^ Throughout this paper, we refer to "selection effects" or "se- 
lection biases" relative to a mass-limited sample, since we arc 
primarily interested in deriving scaling relations as a function of 
cluster mass. 



^ The coverage of REFLEX actually extends slightly into the 
northern hemisphere. In this study, we restrict REFLEX to de- 
clinations < and adjust the sky coverage fraction accordingly. 
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logio ^■"'"'^^^ „ = -13.1931 - 2.31671 logio 2 
ergs-i cm-^ 

-0.10348 (logio (2) 

approximately corresponds to a fixed, intrinsic luminosity 
of 2.5 X 10**/ifo^ ergs"^ in the 0.1-2.4 keV band for our ref- 
erence cosmology. Using this selection criterion, the BCS, 
REFLEX and MACS samples respectively contribute 78, 
126 and 34 clusters to our data set.^ Within the redshift 
and flux range defined above, all three samples are thought 
to be approximately 100 per cent complete and pure.* 

2.2 Follow-up X-ray observations 

Additional information about a subset of the fiux-selected 
clusters is available from follow-up X-ray observations. 
These observations enable more precise measurements of 
cluster luminosities, as well as measurements of morphol- 
ogy, temperature, and gas mass, which is an excellent proxy 
for total mass. 

We incorporate Chandra data for clusters at z > 0.2 
where it is available. Below this redshift, the Chandra Ad- 
vanced CCD Imaging Spectrometer (ACIS) field of view is 
typically too small to allow measurements at rsoo, a canon- 
ical radius at which simulations indicate that the scaling 
relations are well behaved (Evrard et al. 2008). For z < 
0.2, we do, however, incorporate follow-up ROSAT Posi- 
tion Sensitive Proportional Counter (PSPC) observations. 
The ROSAT and Chandra observations were cross-calibrated 
using clusters in the redshift range 0.2 < z < 0.31 that 
were observed with both instruments (Section 2.2.4). De- 
tails of the follow-up observations can be found in Tables 1- 
3; the total good exposure time of the follow-up data, in- 
cluding both ROSAT and Chandra observations, is 3.3 Ms, 
distributed over 94 clusters. The reduction and analysis of 
these data are described in the following sections. 

2.2.1 Chandra data analysis 

The standard level- 1 event files distributed by the Chandra 
X-ray Center (CXC) were reprocessed in accordance with 
CXC recommendations, using the CIAO software package* 

3 Our previous analysis in M08 included 130 REFLEX 
sources using similar selection criteria. Published catalog en- 
tries RX J0117.8-5455 and RX J2251. 7-3206 were eliminated 
from the present study, since high-resolution Chandra imaging re- 
veals their emission to be point-like rather than extended. Source 
RX J0507. 6-0238 satisfies our criteria after correcting its red- 
shift (Dale Koccvski, private communication) from that pub- 
lished in Bohringer et al. (2004), and was therefore added. Three 
other sources (all Abell clusters), RX J0552.8-2103, RX J1959.1- 
3454 and RX J2331. 2-3630, were removed because they do not 
satisfy the flux limit defined in Equations 1 and 2, refiecting 
a small difference between our K-corrcction and that used by 
Bohringer et al. (2004), and the fact that Equation 2 is not pre- 
cisely identical to a fixed luminosity. The BCS and MACS sources 
selected are the same as in MOB. 

* In the sense that the published selection function for each sam- 
ple, the likelihood for cluster detection as a function of redshift 
and flux, is accurate. See additional comments in M08. 

* http://cxc.harvard.edu/ciao/ 



(version 4.1.1, CALDB 4.1.2). This processing includes re- 
moval of bad pixels, corrections for cosmic ray afterglows 
and charge transfer inefficiency, and application of standard 
grade and status filters, using appropriate time-dependent 
gain and calibration products. The extra information avail- 
able for observations taken in VFAINT mode (the majority) 
was used to improve cosmic ray rejection. The data were 
cleaned to remove times of high or unstable background us- 
ing the energy ranges and time bins recommended by the 
CXC. Blank-field data sets made available by the CXC were 
tailored to each observation and cleaned in an identical man- 
ner to the real data. The normalizations of these blank-sky 
files were scaled to match the count rates in the target ob- 
servations measured in the 9.5-12 keV band. 

After identifying and masking point sources, fiat-fielded 
images and background-subtracted surface brightness pro- 
files were prepared in the 0.7-2.0 keV energy band. The 
emissivity in this band is largely insensitive to the gas tem- 
perature, kT, provided kT > 3keV (as is the case for all 
the clusters studied here), making it the preferred energy 
range for determining the gas mass. The center of each clus- 
ter was identified with the centroid of 0.7-2.0 keV emission 
after masking point sources. 

Our spectral analysis is a two-stage process. In the first 
stage, spectra in the 0.8-7.0 keV band were extracted in an 
annulus about the cluster centers. The inner radii of all an- 
nuli were set to 100 kpc in order to prevent relatively cool 
gas in the cores from strongly infiuencing the spectral fits. 
The outer radii were initially chosen to be the radii at which 
the signal-to-noise ratio apparent in the 0.8-7.0 keV surface 
brightness profiles falls to 2. Background spectra were ex- 
tracted from the blank-sky fields for targets at 2; < 0.3, 
and from source-free regions of the detector for z > 0.3 
targets. When required, e.g. due to the presence of strong 
excess soft emission in the field, a model for additional, soft 
thermal emission was included in the spectral modeling of 
the background. Photon-weighted response matrices and ef- 
fective area files were generated for each observation using 
calibration files for the appropriate period. 

The spectral analysis was performed using 
XSPEC.® The spectra were fit to a single-temperature, 
optically thin thermal emission model evaluated 
with the MEKAL code (Kaastra & Mewe 1993; 
Mewe, Kaastra, & Liedahl 1995, incorporating the Fe- 
L calculations of Liedahl, Osterheld, & Goldstein 1995), 
including a model for photoelectric absorption due to 
Galactic hydrogen (Balucinska-Church & McCammon 
1992). Column densities were fixed to the Galactic values 
determined from the Hi survey of Kalberla et al. (2005), 
unless the published value exceeded 10^^ cm^^, in which 
case they were fit as free parameters.'' The abundances 
of all metals were assumed to vary with a common ratio, 

^ http : //heasarc . gsfc.nasa.gov/docs/xanadu/xspec/ 

^ We note that the choice of column density does not have a 
significant impact on our derived luminosities, temperatures and 
masses. In particular, the energy bands used {E > 0.7 keV), and 
the use of locally determined backgrounds for the surface bright- 
ness analysis (Section 2.2.3) and, when possible, for the spectral 
analysis, make our results relatively insensitive to uncertainties 
in the soft X-ray background. Residual uncertainties are within 
the systematic tolerances defined in Section 2.2.3. 
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Table 1. Details of the follow-up observations of BCS clusters. The detector column indicates whether the observation used the ROSAT 
PSPC or the one of the Chandra ACIS CCD cameras. 



Name 


RA {J2000) 


Dec {J2000) 


Date 




Detector 


Mode 


Exposure (ks) 


Abell 2256 


17 


04 02.3 


+78 


38 


14 


1990 


Jun 


17 


PSPC 






16.4 


Abell 1795 


13 


48 53.1 


+26 


35 


38 


1991 
1992 


Jul 
Jan 


01 
04 


PSPC 
PSPC 






15.1 
33.1 


Abell 401 


02 


58 57.2 


+ 13 


34 


46 


1992 
1992 


Jan 
Jul 


23 
30 


PSPC 
PSPC 






5.3 
6.9 


Abell 2029 


15 


10 55.9 


+05 


44 


44 


1992 


Aug 


10 


PSPC 






8.7 


Abell 2255 


17 


12 50.0 


+64 


03 


43 


1993 


Aug 


24 


PSPC 






11.8 


Abell 478 


04 


13 25.2 


+10 


28 


00 


1991 


Aug 


31 


PSPC 






21.9 


Abell 2142 


15 


58 21.1 


+27 


13 


45 


1992 
1992 
1993 


Aug 
Aug 
Jul 


25 
26 
23 


PSPC 
PSPC 
PSPC 






5.9 
4.8 
4.4 


Abell 2244 


17 


02 42.0 


+34 


03 


24 


1992 


Sep 


21 


PSPC 






2.9 


Abell 2034 


15 


10 12.6 


+33 


30 


45 


1993 


Jan 


31 


PSPC 






4.9 


Abell 1068 


10 


40 44.4 


+39 


57 


13 


1992 


Nov 


30 


PSPC 






9.8 


Abell 2204 


16 


32 47.2 


+05 


34 


33 


1992 


Sep 


04 


PSPC 






5.2 


Abell 2218 


16 


35 53.0 


+66 


12 


36 


1991 


May 


25 


PSPC 






31.9 


Abell 1914 


14 


26 01.0 


+37 


49 


37 


1992 


Jul 


20 


PSPC 






6.3 


Abell 665 


08 


30 58.1 


+65 


51 


02 


1991 


Apr 


10 


PSPC 






32.0 


Abell 520 


04 


54 09.0 


+02 


55 


18 


1993 


Sep 


06 


PSPC 






4.8 














2000 


Oct 


10 


ACIS- 


-I 


VFAINT 


9.4 














2003 


Dec 


04 


ACIS- 


-I 


VFAINT 


57.8 














2007 


Jan 


01 


ACIS- 


-I 


VFAINT 


5.1 


Abell 963 


10 


17 03.6 


+39 


02 


52 


1993 


Oct 


28 


PSPC 






9.1 














2000 


Oct 


11 


ACIS- 


-S 


FAINT 


36.3 














2007 


Feb 


18 


ACIS- 


-I 


VFAINT 


5.0 


RX J0439. 0+0520 


04 


39 02.2 


+05 


20 


42 


2000 


Aug 


29 


ACIS- 


-I 


VFAINT 


9.6 














2007 


Nov 


12 


ACIS- 


-I 


VFAINT 


16.2 














2007 


Nov 


15 


ACIS- 


-I 


VFAINT 


7.9 


Abell 1423 


11 


57 17.4 


+33 


36 


40 


2000 


Jul 


07 


ACIS- 


-I 


VFAINT 


9.9 


Zwicky 2701 


09 


52 49.2 


+51 


53 


06 


2001 


Nov 


04 


ACIS- 


-S 


VFAINT 


25.9 


Abell 773 


09 


17 52.7 


+51 


43 


36 


1991 


Oct 


28 


PSPC 






11.0 














2000 


Sep 


05 


ACIS- 


-I 


VFAINT 


10.9 














2003 


Jan 


25 


ACIS- 


-I 


VFAINT 


7.6 














2004 


Jan 


21 


ACIS- 


-I 


VFAINT 


19.8 


Abell 2261 


17 


22 27.0 


+32 


07 


58 


1999 


Dec 


11 


ACIS- 


-I 


VFAINT 


6.7 














2004 


Jan 


14 


ACIS- 


-I 


VFAINT 


23.3 


Abell 1682 


13 


06 50.7 


+46 


33 


30 


2002 


Oct 


19 


ACIS- 


-I 


VFAINT 


1.8 


Abell 1763 


13 


35 19.0 


+40 


59 


59 


1992 


Jun 


23 


PSPC 






12.7 














2003 


Aug 


28 


ACIS- 


-I 


VFAINT 


19.3 


Abell 2219 


16 


40 20.3 


+46 


42 


30 


1993 


Aug 


02 


PSPC 






8.5 














2000 


Mar 


31 


ACIS- 


-S 


FAINT 


42.2 


Abell 2111 


15 


39 41.1 


+34 


25 


07 


1993 


Jul 


23 


PSPC 






6.7 














2000 


Mar 


22 


ACIS- 


-I 


FAINT 


9.5 


Zwicky 5247 


12 


34 22.1 


+09 


47 


05 


2000 


Mar 


23 


ACIS- 


-I 


VFAINT 


9.3 


Abell 267 


01 


52 42.2 


+01 


00 


30 


1999 


Oct 


16 


ACIS- 


-I 


FAINT 


7.1 














2003 


Dec 


07 


ACIS- 


-I 


VFAINT 


17.8 


Abell 2390 


21 


53 37.1 


+ 17 


41 


45 


1993 


Nov 


13 


PSPC 






8.4 














1999 


Nov 


05 


ACIS- 


-S 


FAINT 


8.0 














2000 


Oct 


08 


ACIS- 


-S 


FAINT 


9.8 














2003 


Sep 


11 


ACIS- 


-S 


VFAINT 


74.6 


Zwicky 2089 


no 


UU oD.o 


1 on 
+zU 


oo 


/in 


zUUD 


uec 


23 


Ado- 


T 

-i 


VFAINT 


9.1 


RX J2129. 6+0005 


21 


29 39.7 


+00 


05 


18 


2000 


Oct 


21 


ACIS- 


-I 


VFAINT 


8.7 


RX J0439. 0+0715 


04 


39 00.6 


+07 


16 


00 


1999 


Oct 


16 


ACIS- 


-I 


FAINT 


6.3 














1999 


Oct 


16 


ACIS- 


-I 


VFAINT 


1.6 














2003 


Jan 


04 


ACIS- 


-I 


VFAINT 


19.0 


Abell 1835 


14 


01 01.9 


+02 


52 


40 


1993 
1994 


Jul 
Jul 


03 
05 


PSPC 
PSPC 






6.1 
2.5 














2005 


Dec 


07 


ACIS- 


-I 


VFAINT 


36.3 














2006 


Jul 


24 


ACIS- 


-I 


VFAINT 


39.5 














2006 


Aug 


25 


ACIS- 


-I 


VFAINT 


117.9 


Abell 68 


00 


37 05.9 


+09 


09 


36 


2002 


Sep 


07 


ACIS- 


-I 


VFAINT 


10.0 
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Table 1 - continued 



Name 




RA (J2000) 


Dec (J2000) 


Date 




Detector 


Mode 


Exposure (ks) 


MS J1455. 0+2232 


14 57 15.1 


+22 20 33 


2000 May 19 


ACIS-I 


FAINT 


9.9 










2003 Sep 05 


ACIS-I 


VFAINT 


91.9 










2007 Mar 23 


ACIS-I 


VFAINT 


7.1 


AbcU 697 




08 42 57.6 


+36 21 57 


1991 Apr 10 


PSPC 




9.6 










2002 


Dec 


15 


ACIS-I 


VFAINT 


17.5 


Zwicky 3146 




10 23 39.6 


+04 11 12 


1993 


Nov 


13 


PSPC 




8.1 










2000 May 10 


ACIS-I 


FAINT 


42.4 










2008 


Jan 


18 


ACIS-I 


VFAINT 


31.2 


Abell 781 




09 20 25.3 


+30 30 11 


2000 


Oct 


03 


ACIS-I 


VFAINT 


8.7 


:tails of the follow-up observations of REFLEX clusters (see caption for Table 1). 


Name 




RA (J2000) 


Dec {J2000) 


Date 




Detector 


Mode 


Exposure (ks) 


AT 11 ocrtro 

Abell 35oo 




16 11 58.4 


—31 30 04 


1991 


Jul 


17 


PSPC 




24.0 


A „1 1 or 

Abell oo 




UU 41 oU.i 


— U9 lo Al 


1992 


Jun 


11 


PdPC 




3.2 










1992 


Jul 


01 


PSPC 




8.3 


AK^ill Odd'? 

Abell OD6 1 




on 1 o Q 1 o 


ACi 

-56 49 4/ 


1992 


Oct 


09 


PbPC 




6.3 


AK^ill o^aa 
Abell oz6D 




U4 oi 21.0 


f?-l o^J Ol 

— oi ZD zi 


1992 


Apr 


30 


PoPC 




6.6 










1993 


Aug 


19 


PSPC 




13.2 


AU„11 Q110 

Abell oilz 




r\n 1 coo 
Do I ( DO. 2 


— 44 14 U/ 


1992 


Dec 


17 


PdPC 




7.0 


Au„ll ornv 
Abell 2oy ( 




OQ o c; T n 


1 o nv O'z 
— 12 U ( 2( 


1991 


Nov 


27 


PoPC 




6.7 


Ai 11 onoi 

Abell 3921 




oo Ar\ C "7 o 

22 49 57. (J 


—64 25 53 


1992 


Nov 


15 


PSPC 




11.5 


T\TC Till! O 
Mo Jilil.O 




11 1 A 1 O 

il 14 12. 1 


— 38 11 25 


1993 


Jan 


11 


PdPC 




18.5 


Abell 1689 




11 11 on c 

lo 11 29.0 


m or\ oo 

— Ul 20 28 


1992 


Jul 


18 


PSPC 




13.5 


Abell 2163 




Id 15 4d.U 


f\c oo r ,1 

— (Jd (Jo 54 


1992 


Feb 


28 


PSPC 




3.1 










1992 


Sep 


01 


PSPC 




7.0 










2000 


Jul 


29 


ACIS-I 


VrAlN i 


9.4 










2001 


Jun 


16 


Ado — 1 


VrAllN i 


67.4 


Abell 209 




Ul ol Ov5.1 


— 13 36 48 


2000 


Sep 


09 


Ado — 1 


VrAilN i 


8.7 










2003 


Aug 


03 


ACIS-I 


VrAlN 1 


9.6 


RX J1504.1- 


-0248 


lO U4 U ( .O 


no AQ 1 
~\j2 4o lu 


2004 


Jan 


07 


AOlo — 1 


rAilN 1 


12.0 










2005 


Mar 


20 


ACIS-I 


VrAlN 1 


33.8 


RX J0304.1- 


-3656 


r\'3 r\/1 HQ '3 


oa. t^c on 
— OO OO oU 


2008 


Mar 


16 


AClo — 1 


VrAilN i 


19.9 


RX J0237.4- 


-2630 


r\o or? OT /I 

02 37 27.4 


o '3 r\ o o 

— 26 30 28 


2008 


Mar 


03 


ACIS-I 


VrAlN 1 


7.4 


Abell 2667 




23 51 39.7 


-26 04 60 


1992 


Dec 


05 


PSPC 




2.0 










1992 


Dec 


05 


PSPC 




3.0 










2001 


Jun 


19 


ACIS-S 


VFAINT 


9.4 


RX J0638.7- 


-5358 


06 38 47.3 


-53 58 28 


2008 


Apr 


11 


ACIS-I 


VFAINT 


19.9 


RX J0220.9- 


-3829 


02 20 56.5 


-38 28 52 


2008 


Feb 


29 


ACIS-I 


VFAINT 


19.9 


Abell 521 




04 54 07.4 


-10 13 24 


1999 


Dec 


23 


ACIS-I 


VFAINT 


36.8 










2000 


Oct 


13 


ACIS-S 


VFAINT 


18.6 


RX J0307.0- 


-2840 


03 07 02.0 


-28 39 56 


2008 


Mar 


13 


ACIS-I 


VFAINT 


18.1 


RX J2011.3- 


-5725 


20 11 27.2 


-57 25 10 


2004 


Jun 


08 


ACIS-I 


VFAINT 


22.3 


RX J0232.2- 


-4420 


02 32 17.7 


-44 20 55 


2004 


Jun 


08 


ACIS-I 


VFAINT 


7.9 


RX J0528.9- 


-3927 


05 28 53.3 


-39 28 19 


2004 


Mar 


10 


ACIS-I 


VFAINT 


15.9 


RX J0043.4- 


-2037 


00 43 24.8 


-20 37 24 


2008 


Feb 


02 


ACIS-I 


VFAINT 


19.7 


lES 0657-558 


06 58 27.5 


-55 56 32 


1997 


Fob 


28 


PSPC 




4.5 










2000 


Oct 


16 


ACIS-I 


FAINT 


25.3 










2002 


Jul 


12 


ACIS-I 


VFAINT 


81.1 










2004 


Aug 


10 


ACIS-I 


VFAINT 


21.3 










2004 


Aug 


11 


ACIS-I 


VFAINT 


93.9 










2004 


Aug 


14 


ACIS-I 


VFAINT 


77.5 










2004 


Aug 


15 


ACIS-I 


VFAINT 


31.2 










2004 


Aug 


17 


ACIS-I 


VFAINT 


79.8 










2004 


Aug 


19 


ACIS-I 


VFAINT 


70.7 










2004 


Aug 


23 


ACIS-I 


VFAINT 


22.9 










2004 


Aug 


25 


ACIS-I 


VFAINT 


39.1 


Abell 2537 




23 08 22.0 


-02 11 30 


2004 


Sep 


09 


ACIS-S 


VFAINT 


36.2 
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Table 3. Details of the follow-up observations of Bright MACS clusters (see caption for Table 1). The last two observations listed are 
not used in this paper: the MACS J0358 observation took place after this paper was initially submitted, and we conservatively chose to 
exclude the MACS J2311 observation because it coincides with a long-duration background flare. Their derived properties are nevertheless 
included in Table 6 so that the complete sample of Ebeling et al. (2010) is represented. 



Name RA (J2000) Dec (J2000) Date Detector Mode Exposure (ks) 



MACS 


J2245, 


,0+2637 


22 


45 


04, 


,6 


+26 


38 


04 


2002 


Nov 


24 


ACIS- 


-I 


VFAINT 


13, 


,9 


MACS 


J1131, 


.8-1955 


11 


31 


55, 


,6 


-19 


55 


45 


1993 


Jun 


27 


PSPC 






7 


,2 






















2002 


Jun 


14 


ACIS- 


-I 


VFAINT 


13, 


,1 


MACS 


J0014, 


.3-3022 


00 


14 


18, 


,8 


-30 


23 


18 


1992 


Jun 


16 


PSPC 






13, 


,6 






















2001 


Sep 


03 


ACIS- 


-S 


VFAINT 


23, 


,8 






















2006 


Nov 


08 


ACIS- 


-I 


VFAINT 


12, 


,3 






















2007 


Jun 


10 


ACIS- 


-I 


VFAINT 


37, 


,7 






















2007 


Jun 


14 


ACIS- 


-I 


VFAINT 


23, 


,9 


MACS 


J2140, 


.2-2339 


21 


40 


15, 


,2 


-23 


39 


40 


1993 


Nov 


07 


PSPC 






9, 


.2 






















1999 


Nov 


18 


ACIS- 


-s 


VFAINT 


40, 


.5 






















2003 


Nov 


18 


ACIS- 


-s 


VFAINT 


24, 


.0 


MACS 


J0242, 


,5-2132 


02 


42 


35, 


,9 


-21 


32 


26 


2002 


Feb 


07 


ACIS- 


-I 


VFAINT 


9, 


.6 


MACS 


J1427, 


.6-2521 


14 


27 


39, 


,5 


-25 


21 


03 


2002 


Jun 


29 


ACIS- 


-I 


VFAINT 


14, 


.6 






















2008 


Jun 


11 


ACIS- 


-I 


VFAINT 


26, 


.3 


MACS 


J0547, 


,0-3904 


05 


47 


01, 


,5 


-39 


04 


26 


2002 


Oct 


20 


ACIS- 


-I 


VFAINT 


20, 


.9 


MACS 


J0257, 


,6-2209 


02 


57 


41, 


,3 


-22 


09 


13 


2001 


Nov 


12 


ACIS- 


-I 


VFAINT 


20, 


.5 


MACS 


J2049, 


,9-3217 


20 


49 


55, 


,3 


-32 


16 


49 


2002 


Dec 


08 


ACIS- 


-I 


VFAINT 


22, 


.9 


MACS 


J2229, 


,7-2755 


22 


29 


45, 


,2 


-27 


55 


36 


2002 


Nov 


13 


ACIS- 


-I 


VFAINT 


12, 


.9 






















2007 


Dec 


09 


ACIS- 


-I 


VFAINT 


13, 


.3 


MACS 


J1319, 


,9-1-7003 


13 


20 


07, 


,5 


+70 


04 


37 


2002 


Sep 


15 


ACIS- 


-I 


VFAINT 


8 


.7 


MACS 


J0520, 


,7-1328 


05 


20 


42, 


,2 


-13 


28 


47 


2002 


Feb 


10 


ACIS- 


-I 


VFAINT 


19, 


.0 


MACS 


J1931, 


,8-2634 


19 


31 


49, 


,6 


-26 


34 


34 


2002 


Oct 


20 


ACIS- 


-I 


VFAINT 


13, 


.1 


MACS 


J0035, 


,4-2015 


00 


35 


26, 


,2 


-20 


15 


46 


2003 


Jan 


22 


ACIS- 


-I 


VFAINT 


20, 


.8 


MACS 


J0947, 


,2-1-7623 


09 


47 


13, 


,0 


+76 


23 


14 


2000 


Oct 


20 


ACIS- 


-I 


VFAINT 


11, 


.7 


MACS 


J1115, 


,8-f0129 


11 


15 


51, 


,9 


+01 


29 


55 


2003 


Jan 


23 


ACIS- 


-I 


VFAINT 


8 


.7 






















2008 


Feb 


03 


ACIS- 


-I 


VFAINT 


34, 


.7 


MACS 


J0308, 


,9-1-2645 


03 


08 


56, 


,0 


+26 


45 


35 


2002 


Mar 


10 


ACIS- 


-I 


VFAINT 


23, 


.4 


MACS 


J0404, 


,6-1-1109 


04 


04 


32, 


,7 


+11 


08 


11 


2002 


Feb 


20 


ACIS- 


-I 


VFAINT 


17, 


.6 


MACS 


J1532, 


,8-1-3021 


15 


32 


53, 


,8 


+30 


20 


59 


2001 


Aug 


26 


ACIS- 


-s 


VFAINT 


9 


.4 






















2001 


Sep 


06 


ACIS- 


-I 


VFAINT 


10, 


.0 


MACS 


JOOll, 


,7-1523 


00 


11 


42, 


,8 


-15 


23 


22 


2002 


Nov 


20 


ACIS- 


-I 


VFAINT 


20, 


.8 






















2005 


Jun 


28 


ACIS- 


-I 


VFAINT 


37, 


.2 


MACS 


J0949, 


,8-1-1708 


09 


49 


51, 


,8 


+17 


07 


08 


2002 


Nov 


06 


ACIS- 


-I 


VFAINT 


14, 


.3 


MACS 


J1720, 


,2-1-3536 


17 


20 


16, 


,7 


+35 


36 


23 


2002 


Nov 


03 


ACIS- 


-I 


VFAINT 


19, 


.8 






















2005 


Nov 


22 


ACIS- 


-I 


VFAINT 


29, 


.5 


MACS 


J1731, 


,6+2252 


17 


31 


39, 


,2 


+22 


51 


50 


2002 


Nov 


03 


ACIS- 


-I 


VFAINT 


18, 


.5 


MACS 


J2211, 


,7-0349 


22 


11 


45, 


,9 


-03 


49 


42 


2002 


Oct 


08 


ACIS- 


-I 


VFAINT 


15, 


.2 


TV T A O 

MACS 


J0429 


,6—0253 


04 


29 


36, 


,0 


—02 


53 


06 


2002 


Feb 


07 


ACIS- 


-I 


VFAINT 


22 


.1 


MACS 


J0159, 


,8-0849 


01 


59 


49, 


,4 


-08 


49 


60 


2002 


Oct 


02 


ACIS- 


-I 


VFAINT 


17, 


.4 






















2004 


Dec 


04 


ACIS- 


-I 


VFAINT 


35, 


.3 


MACS 


J2228, 


,5+2036 


22 


28 


32, 


,8 


+20 


37 


15 


2003 


Jan 


22 


ACIS- 


-I 


VFAINT 


19, 


.9 


MACS 


J0152, 


,5-2852 


01 


52 


33, 


,9 


-28 


53 


33 


2002 


Sep 


17 


ACIS- 


-I 


VFAINT 


17, 


.3 


MACS 


J1206, 


,2-0847 


12 


06 


12, 


,3 


-08 


48 


06 


2002 


Dec 


15 


ACIS- 


-I 


VFAINT 


23, 


.5 


MACS 


J0417, 


,5-1154 


04 


17 


34, 


,3 


-11 


54 


27 


2002 


Mar 


10 


ACIS- 


-I 


VFAINT 


11, 


.2 


MACS 


J2243, 


,3-0935 


22 


43 


21, 


,4 


-09 


35 


43 


2002 


Dec 


23 


ACIS- 


-I 


VFAINT 


20, 


.2 


MACS 


J1347, 


,5-1144 


13 


47 


30, 


,8 


-11 


45 


09 


2000 


Mar 


05 


ACIS- 


-s 


VFAINT 


8 


.9 






















2000 


Apr 


29 


ACIS- 


-s 


FAINT 


10, 


.0 






















2001 


May 


10 


ACIS- 


-s 


VFAINT 


89, 


.6 






















2003 


Sep 


03 


ACIS- 


I 


VFAINT 


55, 


.8 


MACS 


J0358, 


,8-2955 


03 


58 


51, 


,2 


-29 


55 


22 


2009 


Oct 


18 


ACIS- 


-I 


VFAINT 


8, 


.7 


MACS 


J2311, 


,5+0338 


23 


11 


35, 


,3 


+03 


38 


25 


2002 


Sep 


07 


ACIS- 


-I 


VFAINT 


6 


.5 



Z, relative to solar values (Anders & Grevesse 1989), 
which was a free parameter in the fit for each cluster. We 
minimized the modified C-statistic (Cash 1979; Arnaud 
1996) rather than , as is appropriate for spectra with few 
counts per bin. Once a best fitting model was identified, 
the parameter space was explored using Markov Chain 



Monte Carlo; these Monte Carlo samples were used to 
propagate the uncertainty due to the fit in the conversion 
of 0.7-2.0 keV count rate to 0.1-2.4 keV flux for each 
observation. We note that, since both the count rate and 
flux in this conversion are at soft energies, the conversion is 
only weakly dependent on temperature and has negligible 
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dependence on metallicity, provided that the temperature 
is sufficiently high {kT > 3keV). 

The second stage of the spectral analysis uses an esti- 
mate for rsoo, obtained as described in Section 2.2.3 below. 
New spectra in the annuli (0.15-0.5)r5oo and (0.15-l)r5oo 
were extracted and fit as above. For observations of nearby 
clusters where a large fraction of the (0.15-l)r5oo region falls 
outside the detector, only the (0.15-0.5)r5oo region was an- 
alyzed; we then estimated the (0.15-l)r5oo temperature us- 
ing the relation To.is-i/To.is-o.s = 0.957 ± 0.009, fit from 
the other clusters. (This result is consistent with that of 
Vikhlinin et al. 2009a, who followed a similar procedure; see 
Section 4.1.1 and Figure 6 of that work.) This final step is 
motivated by the desire to measure a temperature within 
'"500, for consistency with the mass and luminosity measure- 
ments (Section 2.2.3), while still reliably excising the cool 
core, if one is present. We note, however, that this level of 
detail is not entirely necessary; the ~ 4 per cent reduction 
in temperature within rsoo compared with the temperature 
measurement from the first stage of the analysis is within 
statistical errors for most systems and has negligible impact 
on the determination of masses and luminosities. We there- 
fore do not further iterate the determination of rsoo using 
this new temperature estimate. The resulting temperature 
measurements are referred to as kTce in the following, the 
subscript indicating "center-excised", and are listed in Ta- 
bles 4-6. 

We note that the uncertainties on temperatures mea- 
sured within rsoo are necessarily somewhat larger than those 
on our initial estimates from higher signal-to-noise regions. 
Ultimately, however, this has negligible effect on our results, 
since this additional uncertainty is smaller than the system- 
atic allowances associated with instrument calibration and 
mass estimations (Section 2.2.3). 

2.2.2 ROSAT data analysis 

The ROSAT FSVC observations were reduced using the Ex- 
tended Source Analysis Software package of Snowden et al. 
(1994). Time intervals of high particle and scattered so- 
lar X-ray background were excluded, along with periods 
of detector instability. Remaining particle and solar X-ray 
backgrounds were modeled in the nominal energy bands of 
0.7-0.9, 0.9-1.3, and 1.3-2.0 keV (standard ROSAT h&nAs 
R5-R7) and spatially distributed into background maps. 
Exposure maps were generated in each energy band from 
detector maps obtained during the RASS; these exposure 
maps take into account vignetting and detector artifacts. 
Flat-fielded images were created by subtracting the back- 
ground maps from the data and correcting with the expo- 
sure maps. Cluster centroids were determined and surface 
brightness profiles were created, after masking all detectable 
non-cluster sources. Due to the limited energy coverage of 
the PSPC, ROSAT observations are insufficient to deter- 
mine temperatures for the hot clusters under study; instead, 
we adopt temperatures from the literature, measured using 
ASCA (Horner 2001) or Chandra (Vikhlinin et al. 2009a) 
(Tables 4-5). Assuming a typical metallicity, Z = 0.3, we 
generated Monte Carlo samples of the ROSAT count rate 
to fiux conversion factor as above, utilizing response and 
effective area files for each observation generated for the ap- 
propriate region (the central ring) of the PSPC. For the 



high-temperature clusters under study, uncertainties in the 
temperatures and metallicities have a very small impact on 
the measured quantities of interest: luminosity, radius (rsoo) 
and gas mass. Associated systematic uncertainties are well 
within the tolerances marginalized over in Section 2.2.3. 

2.2.3 Gas mass, total mass, and luminosity measurements 

Our analysis uses gas mass as a robust, low-scatter proxy for 
total cluster mass. This choice is motivated by the fact that 
Mgas can be measured robustly, independent of the dynami- 
cal state of the cluster; observations have also indicated con- 
sistency and minimal scatter in the gas mass fraction, /gas = 
Afgas/Mtot, at r25oo for hot, massive clusters (Allen et al. 
2008, hereafter AOS), in agreement with hydrodynamical 
simulations (e.g. Nagai, Vikhlinin, & Kravtsov 2007).* The 
value of /gas at rsoo can be inferred from /gas(r2soo) using 
simulated cluster /gas profiles, accounting for appropriate 
systematic uncertainties, providing a simple conversion from 
Afgas to A'/tot at rsoo. 

Deprojected gas mass profiles were generated for each 
cluster based on its 0.7-2.0 keV surface brightness profile, 
using a modified version of the Cambridge X-ray deprojec- 
tion code of White, Jones, & Forman (1997). Briefiy, the gas 
mass profile is non-parametrically reconstructed using the 
surface brightness profile and a model for the total mass 
distribution, which is adjusted to reproduce the measured 
average temperature. As noted previously, the emissivity 
in the 0.7-2.0 keV band is very weakly dependent on tem- 
perature for the hot clusters in our sample, so the derived 
gas mass is largely insensitive to this procedure; typically, 
changes greater than a factor of two in the predicted tem- 
perature are required to shift the gas mass results by la. 
Uncertainties on the gas mass profiles were determined by 
repeating this procedure for many Monte Carlo realizations 
of the surface brightness profiles. 

Hydrodynamical simulations indicate that systematic 
biases in gas mass measurement due to effects such as as- 
phericity, clumping and projection are small for the general 
cluster population at all radii of interest. Motivated by the 
simulations of Nagai et al. (2007), we incorporate a system- 
atic fractional bias (and associated uncertainty) in gas mass 
measurements at rsoo: B = 0.0325 ± 0.06 (consistent with 
zero). With this systematic allowance, we expect that the 
gas mass values determined by our procedure to be robust. 

The gas masses are converted to total masses using the 
measured gas mass fraction of massive clusters. AOS showed 
that, for hot, massive clusters, the ratio /gas(r2soo) is con- 
sistent with a universal value (see also simulations of, e.g.. 
Eke, Navarro, & Frenk 1998; Nagai, Vikhlinin, & Kravtsov 

* Note that, as detailed in AOS, /gas can only be measured re- 
liably for dynamically relaxed clusters, and is most accurately 
measured at radii r ~ r2soo- For dynamically active systems, 
large systematic variations in measured /gas values are expected 
due to strong, time variable departures from hydrostatic equilib- 
rium (Nagai et al. 2007); the same is true of /gas measurements 
at large radii, regardless of dynamical state. However, the true 
(as opposed to measured) /gas values at r ^ r2soo for all clusters 
are still expected to be similar; that is, mergers effectively disrupt 
hydrostatic equilibrium, but do not affect the gas mass fraction 
at these radii. 
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Table 4. Redshifts and derived properties of BCS clusters from follow-up observations. Temperatures in this work are measured in the 
aperture (0.15— l)r'5oo- Luminosities are in the 0.1— 2.4keV band in the cluster rest frame and have been cross-calibrated to the standard 
of ROSAT (Section 2.2.4). Gas masses measured with Chandra have similarly been calibrated to the iJOSAT standard. Total masses 
are independent of this cross-calibration, but do depend on the value of /gas, which is assumed here to be 0.1104 at r25oo. As discussed 
in the text, recent Chandra calibration updates favor an increase in the measured values of /gas at r2500 by 5-10 per cent, implying 
reductions of 3—5 per cent in rsoo and 7-15 per cent in Mgas,500 a-^d M500, with respect to the values listed here. This dependence is 
accounted for fully in our analysis (see also Paper I). Error bars in this table do not include contributions due to systematic uncertainty 
in /gas, or uncertainty in the overall instrument calibration, whose effects are correlated across all clusters. kT references: [1] Horner 
(2001), [2] Vikhlinin et al. (2009a), [3] this work. 



Name z rsoo A/gas, 500 Afsoo ^500 -^^500,cc 

(Mpc) (IO^^Mq) (lO^^'M©) (10'"'ergs~i) (10** ergs"i) (keV) ref 



Abell 2256 


0, 


.0581 


1, 


.33 


± 


0, 


.08 


0.82 


± 


0. 


,09 


7.2 


± 


1, 


.0 


4.6 


± 


0, 


.4 


3.46 


± 


0.32 


6.90 


± 


0, 


.16 


1 


Abell 1795 


0, 


.0622 


1, 


.22 


± 


0, 


.06 


0.63 


± 


0. 


,07 


5.5 


± 


0, 


.8 


6.0 


± 


0, 


.4 


2.48 


± 


0.16 


6.14 


± 


0, 


.10 


2 


Abell 401 


0, 


.0743 


1, 


.49 


± 


0, 


.08 


1.16 


± 


0. 


15 


10.1 


± 


1, 


.6 


6.8 


± 


0, 


.4 


4.45 


± 


0.29 


7.72 


± 


0, 


.30 


2 


Abell 2029 


0, 


.0779 


1, 


.45 


± 


0, 


.07 


1.07 


± 


0. 


13 


9.3 


± 


1, 


.4 


10.6 




1, 


.0 


4.14 


± 


0.39 


8.22 


± 


0, 


.16 


2 


Abell 2255 


0, 


.0809 


1, 


.24 


± 


0, 


.08 


0.67 


± 


0. 


10 


5.9 


± 


1, 


.0 


2.9 




0, 


.3 


2.39 


± 


0.22 


6.42 


± 


0, 


.16 


1 


Abell 478 


0, 


.0881 


1, 


.48 


± 


0, 


.08 


1.15 


± 


0. 


15 


10.1 


± 


1, 


.6 


13.3 




1, 


.2 


4.88 


± 


0.47 


7.96 


± 


0, 


.27 


2 


Abell 2142 


0, 


.0904 


1, 


.64 


± 


0, 


.08 


1.59 


± 


0. 


19 


13.9 


± 


2, 


.1 


12.4 




0, 


.6 


6.43 


± 


0.35 


10.04 


± 


0, 


.26 


2 


Abell 2244 


0, 


.0989 


1, 


.25 


± 


0, 


.07 


0.70 


± 


0. 


11 


6.2 


± 


1, 


.1 


5.2 


± 


0, 


.5 


2.72 


± 


0.27 


5.37 


± 


0, 


.12 


2 


Abell 2034 


0, 


.113 


1, 


.28 


± 


0, 


.07 


0.77 


± 


0. 


10 


6.7 


± 


1, 


.0 


4.0 


± 


0, 


.4 


2.88 


± 


0.27 


7.15 


± 


0, 


.32 


1 


Abell 1068 


0, 


.1386 


1, 


.04 


± 


0, 


.05 


0.43 


± 


0. 


,06 


3.7 


± 


0, 


.6 


5.6 


± 


0, 


.5 


1.78 


± 


0.18 


3.87 


± 


0, 


.12 


1 


Abell 2204 


0, 


.1511 


1, 


.46 


± 


0, 


.07 


1.18 


± 


0. 


,14 


10.3 


± 


1, 


.5 


17.9 


± 


1, 


.6 


5.30 


± 


0.52 


8.55 


± 


0, 


.58 


2 


Abell 2218 


0, 


.171 


1, 


.28 


± 


0, 


.09 


0.82 


± 


0. 


12 


7.2 


± 


1, 


.2 


5.1 




0, 


.5 


3.35 


± 


0.32 


6.97 


± 


0, 


.37 


1 


Abell 1914 


0, 


.1712 


1, 


.46 


± 


0, 


.10 


1.21 


± 


0. 


14 


10.6 


± 


1, 


.5 


11.3 




1, 


.0 


5.44 


± 


0.55 


9.48 


± 


0, 


.49 


1 


Abell 665 


0, 


.1818 


1, 


.55 


± 


0, 


.07 


1.46 


± 


0, 


16 


12.7 


± 


1, 


.8 


8.6 




0, 


.8 


5.65 


± 


0.52 


8.03 


± 


0, 


.24 


1 


Abell 520 


0, 


.203 


1, 


.50 


± 


0, 


.08 


1.37 


± 


0, 


15 


11.9 


± 


1, 


.6 


8.4 




0, 


.3 


6.41 


± 


0.20 


7.23 


± 


0, 


.23 


3 


Abell 963 


0, 


.206 


1, 


.25 


± 


0, 


.06 


0.78 


± 


0, 


10 


6.8 


± 


1, 


.0 


6.5 


± 


0, 


.2 


3.47 


± 


0.15 


6.08 


± 


0, 


.30 


3 


RX J0439. 0-1-0520 


0, 


.208 


0, 


.92 


± 


0, 


.05 


0.31 


± 


0. 


,04 


2.7 


± 


0, 


.5 


4.4 


± 


0, 


.1 


1.56 


± 


0.06 


4.96 


± 


0, 


.54 


3 


Abell 1423 


0, 


.213 


1, 


.35 


± 


0, 


.10 


1.00 


± 


0. 


,20 


8.7 


± 


2, 


.0 


6.2 


± 


0, 


.4 


3.97 


± 


0.24 


5.75 


± 


0, 


.59 


3 


Zwicky 2701 


0, 


.214 


1, 


.04 


± 


0, 


.06 


0.46 


± 


0. 


,07 


4.0 


± 


0, 


.7 


4.5 


± 


0, 


.2 


1.79 


± 


0.11 


6.75 


± 


0, 


.54 


3 


Abell 773 


0, 


.217 


1, 


.34 


± 


0, 


.06 


0.98 


± 


0. 


,10 


8.6 


± 


1, 


.1 


7.5 


± 


0, 


.2 


4.77 


± 


0.15 


7.37 


± 


0, 


.45 


3 


Abell 2261 


0, 


.224 


1, 


.59 


± 


0, 


.09 


1.65 


± 


0. 


,25 


14.4 


± 


2, 


.6 


12.0 


± 


0, 


.4 


5.63 


± 


0.23 


6.10 


± 


0, 


.32 


3 


Abell 1682 


0, 


.226 


1, 


.50 


± 


0, 


.13 


1.41 


± 


0. 


,34 


12.4 


± 


3, 


.2 


6.5 


± 


0, 


.7 


4.96 


± 


0.61 


7.01 


± 


2, 


.14 


3 


Abell 1763 


0, 


.2279 


1, 


.67 


± 


0, 


.11 


1.94 


± 


0, 


,33 


17.0 


± 


3, 


.4 


10.5 


± 


0, 


.6 


7.29 


± 


0.39 


6.32 


± 


0, 


.40 


3 


Abell 2219 


0, 


.2281 


1, 


.74 


± 


0, 


.08 


2.16 


± 


0, 


,23 


18.9 


± 


2, 


.5 


15.5 


± 


0, 


.8 


9.58 


± 


0.53 


10.90 


± 


0, 


.53 


3 


Zwicky 5247 


0, 


.229 


1, 


.31 


± 


0, 


.14 


0.94 


± 


0, 


18 


8.2 


± 


1, 


.8 


4.3 


± 


0, 


.3 


3.74 


± 


0.24 


5.31 


± 


1, 


.07 


3 


Abell 2111 


0, 


.229 


1, 


.28 


± 


0, 


.11 


0.88 


± 


0. 


,20 


7.8 


± 


1, 


.9 


5.0 


± 


0, 


.3 


3.61 


± 


0.22 


6.51 


± 


0, 


.72 


3 


Abell 267 


0, 


.230 


1, 


.22 


± 


0, 


.07 


0.76 


± 


0. 


,10 


6.6 


± 


1, 


.1 


5.8 


± 


0, 


.2 


3.36 


± 


0.16 


7.13 


± 


0, 


.71 


3 


Abell 2390 


0, 


.2329 


1, 


.61 


± 


0, 


.07 


1.73 


± 


0. 


,17 


15.2 


± 


1, 


.9 


17.3 


± 


0, 


.5 


8.69 


± 


0.29 


10.28 


± 


0, 


.38 


3 


Zwicky 2089 


0, 


.2347 


0, 


.95 


± 


0, 


.04 


0.35 


± 


0, 


,04 


3.1 




0, 


.4 


5.8 




0, 


.3 


1.70 


± 


0.13 


6.55 


± 


1, 


.47 


3 


RX J2129. 6-1-0005 


0, 


.235 


1, 


.28 


± 


0, 


.07 


0.88 


± 


0. 


12 


7.7 


± 


1, 


.2 


9.9 


± 


0, 


.5 


4.17 


± 


0.27 


6.34 


± 


0, 


.62 


3 


RX J0439. 0-1-0715 


0, 


.2443 


1, 


.27 


± 


0, 


.06 


0.85 


± 


0. 


10 


7.4 


± 


1, 


.0 


8.1 


± 


0, 


.3 


4.22 


± 


0.16 


6.59 


± 


0, 


.45 


3 


Abell 1835 


0, 


.2528 


1, 


.49 


± 


0, 


.06 


1.41 


± 


0, 


12 


12.3 


± 


1, 


.4 


21.1 


± 


0, 


.6 


6.73 


± 


0.23 


9.00 


± 


0, 


.25 


3 


Abell 68 


0, 


.2546 


1, 


.27 


± 


0, 


.07 


0.87 




0, 


12 


7.6 




1, 


.2 


6.9 


± 


0, 


.4 


4.42 


± 


0.27 


7.56 


± 


0, 


.97 


3 


MS J1455.0-I-2232 


0, 


.2578 


1, 


.19 


± 


0, 


.06 


0.71 


± 


0, 


,09 


6.2 




1, 


.0 


11.0 




0, 


.3 


3.13 


± 


0.12 


4.54 


± 


0, 


.16 


3 


Abell 697 


0, 


.282 


1, 


.65 


± 


0, 


.09 


1.96 


± 


0, 


,27 


17.1 


± 


2, 


.9 


14.4 




0, 


.8 


8.95 


± 


0.50 


10.93 


± 


1, 


.11 


3 


Zwicky 3146 


0, 


.2906 


1, 


.35 


± 


0, 


.06 


1.08 


± 


0, 


10 


9.4 




1, 


.2 


19.1 


± 


0, 


.7 


5.82 


± 


0.26 


8.38 


± 


0, 


.44 


3 


Abell 781 


0, 


.2984 


1, 


.26 


± 


0, 


.07 


0.90 


± 


0, 


14 


7.9 


± 


1, 


.4 


6.0 


± 


0, 


.4 


5.10 


± 


0.32 


7.55 


± 


1, 


.03 


3 



2007). Given a model for the gas mass fraction profile, the 
gas mass profiles can be used to determine total masses. 
Ma, and associated radii, ta, corresponding to mean en- 
closed overdensity A, by solving the implicit equation 



M(rA) = 



_ 47rApcr(2:)ri 

(l + B)/gas(rA) " 3(1 + B) ' 



(3) 



where per (2) is the critical density at the redshift of the 
cluster, for va and M(rA). Adopting the value /gas (f 2500) = 
0.1104 from AOS, we verified that our analysis produces esti- 
mates of r25oo and A/2500 that are compatible with the more 
detailed results of that work for the systems in common. 

To obtain mass estimates at larger radii (A — 500) 
from the gas mass profiles, we need an estimate of /gas at 



this radius. The increased prevalence of departures from hy- 
drostatic equilibrium at these larger radii (e.g. Nagai et al. 
2007) prevent us from measuring /gas (''500) directly, even 
for dynamically relaxed systems. Instead, we use the predic- 
tion of simulations (e.g. Eke et al. 1998; Grain et al. 2007) 
that hot, massive clusters should display approximately self- 
similar /gas profiles at radii r ^ r25oo (approximately one 
quarter of the virial radius), independent of their dynami- 
cal state. Using the simulations of Eke et al. (1998), we fit 
a power-law model to fgasir) between r25oo and rsoo, find- 
ing a logarithmic slope of 0.048 ± 0.030. Assuming the total 
mass follows a Navarro, Frenk, & White (1997) profile with 
a typical concentration parameter c = 4 (Zhao et al. 2003, 
2009; Gao et al. 2008), the ratio r5oo/r25oo is approximately 
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Cluster growth: X-ray scaling relations 

Table 5. Rcdshifts and derived properties of REFLEX elusters from follow-up observations (see caption for Table 4). 
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Name Z rsOO A/gas, 500 Afsoo ^500 -^^500, CO kTcc kT 

(Mpc) (IO^^Mq) (lO^^'Mo) (lO'^'ergs-i) (10** ergs^^) (keV) ref 



Abell 3558 




0, 


.048 


1, 


.28 


± 


0.10 


0.73 


± 


0, 


.09 


6.4 


± 


1, 


.0 


3, 


.7 




0.3 


2, 


.60 


± 


0.24 


5.51 


± 


0, 


.10 


1 


Abell 85 




0, 


.0557 


1, 


.33 


± 


0.06 


0.82 


± 


0, 


.09 


7.2 


± 


1, 


.0 


5, 


.7 




0.4 


3, 


.04 


± 


0.20 


6.45 


± 


0, 


.10 


2 


Abell 3667 




0, 


.0557 


1, 


.57 


± 


0.10 


1.35 


± 


0, 


.21 


11.8 




2, 


.2 


5, 


.8 


± 


0.5 


4, 


.38 


± 


0.40 


6.33 


± 


0, 


.06 


2 


Abell 3266 




0, 


.0602 


1, 


.45 


± 


0.07 


1.06 




0, 


.13 


9.2 


± 


1, 


.4 


4, 


.9 


± 


0.3 


3, 


.77 


± 


0.24 


8.63 


± 


0, 


.18 


2 


Abell 3112 




0, 


.0752 


1, 


.10 


± 


0.05 


0.47 


± 


0, 


.06 


4.1 


± 


0, 


.6 


4, 


.4 


± 


0.4 


1, 


.73 


± 


0.16 


4.28 


± 


0, 


.09 


1 


Abell 2597 




0, 


.0852 


0, 


.97 


± 


0.06 


0.33 


± 


0, 


.05 


2.9 


± 


0, 


.5 


4, 


.3 


± 


0.4 


1, 


.29 


± 


0.13 


3.58 


± 


0, 


.07 


1 


Abell 3921 




0, 


.094 


1, 


.20 


± 


0.09 


0.62 


± 


0, 


.09 


5.4 


± 


0, 


.9 


3, 


.1 


± 


0.3 


2, 


.23 


± 


0.21 


5.07 


± 


0, 


.17 


1 


MS J1111.8 




0, 


.1306 


1, 


.12 


± 


0.09 


0.52 


± 


0. 


.07 


4.5 


± 


0, 


.7 


3, 


.2 




0.3 


2, 


.29 


± 


0.22 


5.79 


± 


0, 


.22 


1 


Abell 1689 




0, 


.1832 


1, 


.45 


± 


0.07 


1.21 


± 


0, 


.14 


10.5 


± 


1, 


.5 


13, 


.6 


± 


1.2 


5, 


.72 


± 


0.57 


9.15 


± 


0, 


.35 


1 


Abell 2163 




0, 


.203 


2, 


.22 


± 


0.10 


4.40 


± 


0, 


.45 


38.5 


± 


5, 


.0 


28, 


.7 


± 


1.1 


17, 


.77 


± 


0.67 


12.27 


± 


0, 


.90 


3 


Abell 209 




0, 


.206 


1, 


.53 


± 


0.08 


1.44 


± 


0, 


.18 


12.6 


± 


1, 


.9 


8, 


.6 


± 


0.3 


5, 


.80 


± 


0.22 


8.23 


± 


0, 


.66 


3 


RX J1504.1- 


-0248 


0, 


.2153 


1, 


.46 


± 


0.06 


1.25 


± 


0, 


.13 


11.0 


± 


1, 


.4 


27, 


.6 




1.0 


5, 


.74 


± 


0.26 


8.00 


± 


0, 


.44 


3 


RX J0304.1- 


-3656 


0, 


.2192 


1, 


.07 


± 


0.07 


0.51 


± 


0, 


.08 


4.4 


± 


0, 


.9 


3, 


.0 


± 


0.2 


2, 


.05 


± 


0.12 


6.27 


± 


0, 


.76 


3 


RX J0237.4- 


-2630 


0, 


.2216 


1, 


.16 


± 


0.08 


0.65 


± 


0, 


.11 


5.6 


± 


1, 


.1 


5, 


.6 


± 


0.3 


2, 


.59 


± 


0.18 


6.65 


± 


1, 


.26 


3 


Abell 2667 




0, 


.2264 


1, 


.36 


± 


0.07 


1.03 


± 


0. 


.14 


9.0 


± 


1, 


.5 


12, 


.6 


± 


0.6 


4, 


.97 


± 


0.31 


7.58 


± 


0, 


.73 


3 


RX J0638.7- 


-5358 


0, 


.2266 


1, 


.42 


± 


0.06 


1.18 


± 


0, 


.13 


10.3 


± 


1, 


.4 


11. 


.2 


± 


0.6 


5, 


.48 


± 


0.31 


9.45 


± 


0, 


.95 


3 


RX J0220.9- 


-3829 


0, 


.228 


1, 


.07 


± 


0.08 


0.50 


± 


0, 


.11 


4.4 


± 


1, 


.1 


4. 


.0 


± 


0.2 


2, 


.09 


± 


0.13 


5.17 


± 


0, 


.53 


3 


Abell 521 




0, 


.2475 


1, 


.46 


± 


0.07 


1.31 


± 


0. 


.16 


11.4 


± 


1. 


.7 


7. 


.1 


± 


0.3 


5, 


.80 


± 


0.21 


6.21 


± 


0, 


.28 


3 


RX J0307.0- 


-2840 


0, 


.2537 


1, 


.28 


± 


0.07 


0.90 




0. 


.13 


7.8 


± 


1. 


.4 


7. 


.3 


± 


0.4 


3, 


.73 


± 


0.22 


9.60 


± 


1, 


.49 


3 


RX J2011.3- 


-5725 


0, 


.2786 


0. 


.95 


± 


0.07 


0.37 


± 


0. 


.07 


3.3 


± 


0, 


.7 


3. 


.8 


± 


0.2 


1, 


.67 


± 


0.12 


3.23 


± 


0, 


.34 


3 


RX J0232.2- 


-4420 


0, 


.2836 


1, 


.49 


± 


0.10 


1.45 


± 


0, 


.25 


12.7 


± 


2, 


.5 


13, 


.3 


± 


0.7 


6, 


.57 




0.44 


10.06 


± 


2, 


.31 


3 


RX J0528.9- 


-3927 


0, 


.2839 


1, 


.52 


± 


0.07 


1.52 


± 


0, 


.17 


13.3 


± 


1, 


.8 


11, 


.6 


± 


0.6 


6, 


.92 


± 


0.39 


7.80 


± 


0, 


.85 


3 


RX J0043.4- 


-2037 


0, 


.2924 


1, 


.28 


± 


0.07 


0.92 


± 


0. 


.14 


8.1 


± 


1. 


.4 


7, 


.7 


± 


0.4 


4, 


.75 


± 


0.29 


7.59 


± 


0. 


.77 


3 


lES 0657-558 


0, 


.2965 


1, 


.81 


± 


0.07 


2.61 


± 


0. 


.24 


22.8 


± 


2. 


.8 


21, 


.7 


± 


0.3 


13, 


.61 


± 


0.26 


11.70 


± 


0. 


.22 


3 


Abell 2537 




0, 


.2966 


1, 


.23 


± 


0.06 


0.82 


± 


0, 


.11 


7.2 


± 


1. 


.1 


6, 


.4 


± 


0.4 


3, 


.87 


± 


0.29 


7.63 


± 


0. 


.86 


3 



2.3, yielding a scaled, mean gas mass fraction estimate of 
/gaa(r5oo) ~ 0.115 from the AOS value of /gas(r250o)- We 
note that this procedure is insensitive to the assumed con- 
centration parameter, since (r5oo/''250o)'^ °''* is a weak func- 
tion of c, varying at only the per cent level between c = 2 
and c = 10. The systematic uncertainty associated with the 
concentration is thus small in comparison with the other 
allowances described below, as well as the statistical uncer- 
tainties. 

The dominant systematic uncertainty affecting our to- 
tal mass measurements is associated with the measurement 
/gas (f 2500 )■ For the purposes of obtaining mass estimates at 
''500, we adopted /gas(''250o) = 0.1104 as a reference value, 
similar to the cosmological reference values h = 0.7 and 
f2m = 0.3. The systematic uncertainty on the universal /gas 
value is fully accounted for in the analysis by incorporat- 
ing and marginalizing over the detailed model of AOS (in- 
cluding systematic allowances at the 10-20 per cent level; 
see Section 2.3) and modeling the dependence of measured 
quantities on /gas in a similar manner to the cosmological 
reference parameters (see Paper I). Systematic, cluster-to- 
cluster scatter in /gas is undetected in the sample of large, 
dynamically relaxed clusters studied in AOS; we nevertheless 
incorporated an additional systematic uncertainty in the in- 
dividual mass estimates at the 5 per cent level to account 
both for possible intrinsic scatter and residual calibration or 
modeling uncertainties. This uncertainty is reflected in the 
error bars listed in Tables 4-6. 

Luminosities in the RASS band (0.1-2.4 keV) within 
^500 were calculated by integrating the surface brightness 
profiles, using the count rate to flux conversions computed 
from the spectral analysis. This integration was accom- 



plished by fitting the profile outskirts with a power-law 
cluster emission plus constant background model; the lumi- 
nosities were background subtracted and, where necessary, 
the power-law model was used to extrapolate into the low 
signal-to-noise regime. Uncertainties from the profile mod- 
eling were propagated by Monte Carlo sampling parameters 
of the fits, and were combined with a 5 per cent systematic 
allowance accounting for residual observation-to-observation 
calibration uncertainties. The mass and luminosity measure- 
ments appear in Tables 4-6. 

We also calculated projected luminosities for each clus- 
ter within the annulus 0.15 < r/rsoo < 1. The procedure is 
the same as above, except that the integration of the surface 
brightness profile begins at r = 0.15r5oo rather than r = 0. 
We note that, for the z < 0.2 clusters where we use ROSAT 
data, the resolution of the PSPC is sufficient to perform this 
center-excision. The resulting luminosities, denoted L5oo,cc, 
are investigated in Section 6, but do not contribute to the 
main results of this paper, or the cosmological constraints 
in Paper I. 



2.2.4 Instrument cross-calibration 

There are 11 clusters in our sample for which both ROSAT 
and Chandra follow-up observations were analyzed (see Ta- 
bles 1-3). These clusters were used to cross-calibrate the two 
instruments, although in the end only the Chandra data for 
these 11 clusters were used to determine the properties listed 
in Tables 4-6. For the purposes of the cross-calibration, the 
masses and luminosities from ROSAT were estimated using 
centroids, temperatures and metallicities determined from 
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Table 6. Redshifts and derived properties of Bright MACS clusters from follow-up observations (see caption for Table 4). The last two 
entries represent clusters from Ebeling et al. (2010) whose follow-up observations were not used in this paper, but which we include in 
this table for completeness; for these clusters we did not calculate the center-excised luminosity. 



Name Z rsoo Mgaa,500 ^^500 -^-SOO -^'SOO.cc 

(Mpc) (IQI^Mq) (10" Mq) (10*-*ergs-i) (lO^^ergs"!) (keV) ref 



MACS 


J2245, 


.04 


-2637 


0.301 


1, 


.08 


± 


0, 


.07 


0, 


.56 


± 


0, 


,10 


4, 


,9 


± 


1, 


.0 


7, 


.6 


± 


0, 


,4 


3, 


,09 


± 


0.24 


5, 


.47 


± 


0, 


.58 


3 


MACS 


J1131, 


.8- 


-1955 


0.306 


1, 


.69 


± 


0, 


.09 


2, 


.16 


± 


0, 


.28 


18, 


,9 


± 


2, 


.9 


13, 


.1 


± 


0, 


,7 


8, 


,48 


± 


0.49 


9, 


.35 


± 


1, 


.67 


3 


MACS 


J0014, 


.3- 


-3022 


0.308 


1, 


.65 


± 


0, 


.07 


2, 


.01 
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Chandra and the comparisons were made at a fixed radius, 
the C/iandra-determined value of rsoo- 

Using CIAO version 4.1.1 and CALDB 4.1.2, the lumi- 
nosities determined from Chandra data are systematically 
larger than those from ROSAT by 14 per cent, with a scat- 
ter of 7 per cent. This disagreement is predominantly a con- 
sequence of the 21 January 2009 Chandra calibration up- 
date, which significantly changed the ACIS effective area 
at soft energies. A previous analysis using CALDB version 
3.4.3 and CIAO 3.4 (which also did not incorporate correc- 
tions for dead area in the front-illuminated chips caused by 
cosmic rays, now accounted for by default) produced near 
perfect agreement between Chandra and ROSAT luminosi- 
ties. Since our ultimate goal is to relate cluster masses to 
T^OSAr survey fluxes, we apply this 14 per cent correction 
to the Chandra luminosities in Tables 4-6. Our results on the 
luminosity-mass and temperature-luminosity relations are 
thus for i?05j4r-calibrated luminosities; (7/iancira-calibrated 
relations can be derived simply by adjusting the normaliza- 
tions of the relations (Sections 3 and 5) appropriately. 

Our Chandra gas mass estimates, which to good ap- 
proximation scale as the square root of the luminosity, were 



correspondingly larger by ~ 7 per cent, with a scatter of 
2 per cent. We again calibrated the Chandra gas masses to 
the i?05'j4r standard, which agrees extremely well with pre- 
vious versions of the Chandra calibration. In this case, the 
correction is motivated by the need for consistency between 
gas masses measured in this study with those of AOS, which 
predate the change to the ACIS effective area. In taking the 
ratio /gas/Afgaa (Equation 3), with /gas from AOS and Mgaa 
measured in this study, any common systematic bias in gas 
mass cancels, leaving our estimates of the total masses un- 
biased. We note that the 21 January 2009 Chandra calibra- 
tion update additionally affects temperature measurements, 
which is relevant for the AOS determination of /gas; we ad- 
dress this issue in Section 2.3. 

2.3 /gas data 

Since our analysis uses gas mass as a proxy for total cluster 
mass, we must simultaneously constrain both the gas mass 
fraction and the growth of structure. For that reason, our 
standard analysis includes the /gas analysis of AOS; however, 
we use only the 6 lowest-redshift clusters from that paper 
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{z < 0.15). Those 6 clusters are sufficient to constrain the 
gas mass fraction at low redshift, without themselves pro- 
ducing a direct constraint on the expansion of the Universe 
(and thus unduly influencing our cosmological results). This 
simultaneous analysis incorporates the full /gas model de- 
tailed in AOS, and includes generous systematic allowances 
for instrument calibration (10 per cent), non-thermal pres- 
sure support (10 per cent, Nagai et al. 2007), the depletion 
of baryons in clusters with respect to the cosmic mean (20 
per cent), and evolution with redshift of the baryonic and 
stellar content of clusters (10 and 20 per cent). 

The raw data analysis in AOS predates a recent (21 Jan- 
uary 2009) update to the Chandra ACIS effective area at 
soft energies. As discussed in Section 2.2.4, the effect of the 
calibration on the AOS gas masses can be compensated by 
calibrating our own gas masses consistently. However, the 
calibration update has an additional effect on the determi- 
nation of /gas; the change in measured temperatures due to 
the relative effective area at low and high energies results in 
a corresponding change in total mass estimates. In our tests, 
this effect results in an increase in /gas values for massive 
clusters at the 5-10 per cent level; we have accounted for 
this correction by shifting the center of the "Chandra cal- 
ibration" nuisance parameter in the AOS model such that 
the preferred value of /gas increases correspondingly by 10 
per cent. The Gaussian prior on this parameter has itself 
a width of 10 per cent, so the systematic allowance encom- 
passes both the old value and the new value expected from 
the calibration update. 

2.4 Other data 

As discussed above, our analysis of the XLF data in- 
cludes 6 low redshift {z < 0.15) /gas clusters from AOS. 
As in MOS, when analyzing the X-ray cluster data alone, 
we also incorporate Gaussian priors on the Hubble con- 
stant, h — 0.72 ± O.OS, based on the Hubble Key project 
(Freedman et al. 2001), and on the mean baryon density, 
Q^yh^ = 0.0214 ± 0.002, based on big bang nucleosynthesis 
studies (Kirkman et al. 2003). These priors are not needed 
when cosmic microwave background (CMB) data are also 
included in the analysis; see below. 

In Section 5, we present results from a combination 
of XLF and other independent cosmological data, includ- 
ing cluster gas mass fractions (AOS, including the full 42 
cluster data set), 5-year CMB data from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP5, Dunkley et al. 2009), 
the type la supernova (SNIa) analysis of Kowalski et al. 
(200S, including their treatment of systematic errors), 
and the Baryon Acoustic Oscillation (BAO) analysis of 
Percival et al. (2007). More details on these data and their 
use are reviewed in Paper I. 

When fitting CMB data, we allow the scalar spectral in- 
dex. Us, and optical depth to reionization, r, to vary as free 
parameters, and marginalize over a plausible range in the 
amplitude of the Sunyaev-Zel'dovich signal due to galaxy 
clusters (0 < Asz < 2; introduced by Spergel et al. 2007). 
The combination of CMB and /gas data places tight con- 
straints on both h and fly^h^ in addition to other parame- 
ters of interest (see AOS), so the Hubble Key project and big 
bang nucleosynthesis priors are unnecessary in analyses of 
the combined data sets. This applies as well to the combi- 



nation of CMB and XLF data, since we always use a low-z 
subset of the /gas data to calibrate the cluster mass scale for 
the XLF analysis. 



3 SCALING RELATION MODEL 

The model used for the joint cosmology and scaling rela- 
tion work is described in full in Paper I. Here we review the 
basic scaling relation model and describe various possible 
extensions to it that are investigated in Section 5.3. 

In Paper I, the nominal values of luminosity and tem- 
perature for a given mass are modeled as power laws: 

(iim)) = /3r+/3rm, (4) 
(t{m)) = /?r + /3rm. 

Here the terms 

(} _ ] I -^500 \ 

^ " l^i5(«)10«ergs-i ) ' 

include factors of the normalized Hubble parameter, E{z) = 
H{z)/Ho, expected for "self-similar" evolution of the rela- 
tions with redshift. These factors appear because the rela- 
tion is defined in terms of quantities measured within ra- 
dius rsoo, and follow from the dependence of r^oo on the 
evolving critical density of the Universe (e.g. Kaiser 19S6; 
Bryan & Norman 199S).® The intrinsic scatter in £ and t at 
fixed m is non-negligible, and can be modeled most simply 
as bivariate normal, described by marginal scatters, atrn and 
atm, and a scatter correlation coefficient, pum- We refer to 
this seven parameter model (/3o'", /^i'", ^cT, Pi"^, (^im, c^tm, 
pitm) as the "minimal" scaling relation model. 

This model can be naturally extended in several ways. 
One possibility is to allow departures from self-similar evo- 
lution in the nominal relations. We write 

{l{m)) = /^r-f /jrm + /?^logio(C), (6) 
(t(m)> = /?r-h/3rm + /?nog,o(C), 

where is a function of redshift, commonly chosen to be 
1-1- 2: (e.g. Ettori et al. 2004; Morandi, Ettori, & Moscardini 
2007; MOS) or E{z) (e.g. Maughan 2007; Vikhhnin et al. 
2009a). 

Wc similarly consider the possibility of evolution in the 
intrinsic scatter, making each scatter parameter {aim, cFtm 
and pitm) a function of redshift of the form 

(yim{z) = <7lm[l + a'lm{C, - l)\ ■ (7) 

A third potentially interesting avenue is to allow asym- 
metry in the intrinsic scatter. We consider generalizing the 
bivariate normal distribution used to describe the scatter 
to the bivariate skew-normal distribution (for reference, see 

^ We note that evolution of this form is a generic prediction of 
gravitationally driven collapse, appearing ubiquitously in simula- 
tions (e.g. Evrard et al. 2008), and docs not require the assump- 
tion of virial equilibrium. 
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Azzalini 1985; Azzalini & Valle 1996; Azzalini & Capitanio 
1999). For simplicity, we only consider skewness in the 
marginal luminosity-mass relation, parametrized by the 
shape parameter \im- 



4 ANALYSIS 

The analysis method used in this paper is described in 
full in Paper I. Uniquely, it provides a rigorous and fully 
self-consistent means of producing simultaneous constraints 
on both the cluster scaling relations and cosmological pa- 
rameters; by this we mean that the likelihood of the full 
data set (survey + follow-up observations) is derived from 
first principles, ensuring that the covariance among all the 
model parameters is fully captured and that selection bi- 
ases^" (Malmquist bias and Eddington bias) are properly 
accounted for. The results are determined by sampling of 
the posterior distribution using Markov Chain Monte Carlo, 
marginalizing over priors encoding all appropriate system- 
atic uncertainties (see Paper I). The most important aspects 
in this context are that (1) the sample selection functions^^ 
are included; (2) the cluster mass function is taken into ac- 
count; and (3) the correlation of measurement errors is prop- 
erly incorporated into the likelihood. As reviewed in Ap- 
pendix A, the first two points above are crucial for a proper 
accounting of Malmquist and Eddington biases. The above 
factors, as well as the large size and high completeness and 
purity of the cluster samples, and the quality of the follow-up 
data, lead to significant improvements over previous work. 

Consistent with the method described above, all 238 
clusters meeting our selection criteria are used in the anal- 
ysis presented here. However, the 94 clusters with follow-up 
data naturally have more influence on the scaling relation 
results (for more details, see Paper I, Section 4.1.3). 

In this paper, the only cosmological model we con- 
sider is the simple, spatially flat, cosmological constant 
model (ACDM) which is known to provide a good fit to 
CMB and kinematic data (e.g. AOS; Kowalski et al. 2008; 
Dunkley et al. 2009). Our results for the cluster scaling re- 
lations are marginalized over the parameters of this cosmo- 



Sincc wc arc concerned with selection biases relative to a mass- 
selected sample, the intrinsic dispersion between X-ray luminosity 
and mass is an important (indeed, the dominant) contributing 
factor, in addition to survey flux measurement errors (see also 
Appendix A). 

By "selection function" , we refer to the probability that a clus- 
ter be detected by a survey and included in the data set, as a 
function of that cluster's physical properties (redshift, flux etc.). 
See Section 4.1 of Paper I. 

While the mass function clearly must be taken into account 
to properly handle Eddington and Malmquist biases (see Ap- 
pendix A), there is a second, independent effect for which it is 
relevant. Generically, but especially in the presence of intrinsic 
scatter, the distribution of covariates (in this case m) is signifi- 
cant. For example, in methods based on least squares, the data at 
the ends of the range of covariates are disproportionately influen- 
tial in determining the slope; thus, even in schemes that are osten- 
sibly modified to include intrinsic scatter (e.g. Press et al. 1992), 
the assumption of uniformly distributed covariates can produce 
biased results (see Kelly 2007 and references therein). 



Table 7. Marginalized 68.3 per cent confidence intervals on the 
parameters of the minimal (7 parameter) scaling relation model 
from XLF data alone, and from the combination of all the cosmo- 
logical data sets (Section 2.4). Tildes indicate parameters describ- 
ing the bolometric luminosity-mass (as opposed to 0.1-2.4 keV 
luminosity) relation; the bolometric results are derived from the 
0.1-2.4 keV results using the temperature-mass relation, indi- 
cates temperature measured in the full r < rsoo region of each 
cluster (without excising the center). 
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crtm 


0.056 


± 


0. 


,008 


0.055 


±0, 


.008 


scatter correlation 


Pltm 


0.10 


± 


0. 


19 


0.09 


±0, 


.19 


^bol^™ normalization 


Po 


1.26 


± 


0. 


17 


1.23 


±0, 


.12 


£bol— m. slope 




1.59 


± 


0. 


,09 


1.63 


±0, 


.06 


tci-m normalization 


otm 
Po.ci 


0.90 


± 


0, 


,04 


0.91 


±0, 


.03 


tci-m slope 


otm 
Pl,ci 


0.46 


± 


0, 


,04 


0.46 


±0, 


.04 


tci-m marginal scatter 




0.059 


± 


0. 


,007 


0.06C 


l±C 


1.07 



logical model, which are enumerated in Paper I. Constraints 
on cosmology from our data are also presented in that work. 



5 RESULTS ON THE SCALING RELATIONS 
5.1 Constraints on the minimal model 

5.1.1 The joint luminosity-temperature-mass relations 

One-dimensional, marginalized constraints on the parame- 
ters of the minimal scaling relation model from the XLF 
data are presented in Table 7. Joint constraints on the 
most degenerate parameter pairs are shown in Figure 2; 
the only significant degeneracy for this simple model is be- 
tween the normalizations of the nominal luminosity-mass 
and temperature-mass relations, both of which strongly cor- 
relate with /gas and cosmological parameters. ^'^ 

The luminosity-mass and temperature-mass data and 
the best-fitting nominal relations from the analysis are 
shown in Figure 3. For the purposes of this visual com- 
parison, the data have been adjusted by the E{z) factors 
appropriate to the best-fitting cosmology and to account 
for the cosmological dependence of rsoo (a smaller correc- 
tion; see Section 3.4 of Paper I). For the luminosity-mass 
relation in the 0. 1-2.4 keV band, we measure normalization 
/3o™ = 0.82 ± 0.11, slope ^f™ = 1.29 ± 0.07, and scatter 
aim = 0.180 ± 0.019 (~ 40 per cent). 

The apparent ofi'set in normalization of the best fit (red 
line) from the luminosity-mass data refiects the dramatic 
effect of Malmquist and Eddington biases on the data set. 

The normalization and slope of each nominal relation arc not 
strongly degenerate because we have balanced the data around 
an appropriate pivot (10^^ A/q; Equation 5). 
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Figure 2. Joint 68.3 and 95.4 per cent confidence regions on parameters of the minimal (7 parameter) scaling relation model from XLF 
data (purple) and from the combination with independent cosmological data (gold, Section 2.4). Left: normalizations of the nominal 
luminosity— mass and temperature-mass relations. The degeneracy is principally due to the fact that both parameters correlate with 
the overall mass scale set by /gas- Center: slope and marginal intrinsic scatter of the luminosity-mass relation. Right: same for the 
temperature-mass relation. The combination with other data is not shown in this panel, since the results are very similar (Table 7). 





Figure 3. Mass, luminosity and temperature data for our cluster sample. Note that the data are adjusted with factors of E{z) appropriate 
for the best-fitting cosmology (Paper I), not the reference cosmology. Also shown are the best- fitting, nominal luminosity— mass and 
temperature— mass power-law relations (solid, red lines), and the results of fits to the data with the BCES Y\X and X\Y methods (purple 
lines, respectively dashed and dot-dashed). The best-fitting values of the intrinsic scatter in each relation, t7£„i and crtm, are indicated 
by the orange points with error bars in the lower-right corner of each panel. The seemingly low normalization of our luminosity— mass fit 
compared with the data reflects the effects of Malmquist and Eddington biases on the sample, which are not accounted for in the BCES 
fits. Plotted error bars do not include the contributions due to systematic uncertainty in /gas, or uncertainty in the overall instrument 
calibration, although these are included in the analysis. The correlation between luminosities and masses measured from the same X-ray 
data are accounted for (see Paper I); this small correlation (p fa 0.3) is illustrated by the blue ellipse in the left panel, which shows the 
(joint, 68.3 per cent confidence) error ellipse for one cluster. There is no correlation between measured temperatures and masses, using 
our procedure (Section 2.2). 



Specifically, the flux limits of our cluster samples are liigli 
enougii tiiat the data set consists largely of sources that 
are of average luminosity or greater, given their masses. 
Pacaud et al. (2007) observed a similar effect in an analysis 
of XMM-Newton Large Scale Structure (XMM-LSS) clus- 
ters that accounted for the survey selection function, as did 
Stanek et al. (2006), who incorporated both the selection 
function and cluster mass function into an analysis of the 
REFLEX sample. These results underscore the fact that se- 



lection effects must be fully accounted for in order to accu- 
rately recover the scaling relations from cluster data. 

In contrast, naively applying the commonly used BCES 
Y\X or X\Y methods of Akritas & Bershady (1996) to the 
plotted data results in the purple dashed and dot-dashed 
lines in Figure 3, respectively. For the luminosity-mass rela- 
tion (left panel), both BCES fits, which do not account for 
selection bias, have slopes and normalizations that are incon- 
sistent with the true best fit (red, solid line). The discrep- 
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ancy is significant for cosmological applications, since the 
luminosity-mass relation is influential in determining the to- 
tal number of cluster detections, and is thus degenerate with 
and its. We have similarly flt the data in Figure 3 with 
other simple regression methods based on modified least- 
squares techniques (the fitexy routine of Press et al. 1992) 
and Bayes theorem^** (Kelly 2007), and find that they pro- 
duce results very similar to the BCES Y\X and X\Y meth- 
ods, depending on which variable is assigned to be the co- 
variate. 

For the temperature-mass relation (right panel of Fig- 
ure 3; using center-excised temperatures, kTcc), we find nor- 
malization Pl'^ = 0.88 ± 0.05, slope /J^™ = 0.48 ± 0.04, and 
scatter atm ~ 0.056 ± 0.008 (~ 13 per cent scatter). Since 
this relation is not as influenced by selection bias as the 
luminosity-mass relation is, the BCES flts are much closer 
to our best fit in normalization in this case. 

The slope /Sf™ can be simply converted to a bolomet- 
ric luminosity-mass slope, /Sf™, using the temperature-mass 
relation. (We use a tilde to distinguish parameters describ- 
ing the bolometric luminosity relation.) The low scatter 
of the temperature-mass relation, together with the fact 
that the ratio of 0.1-2.4 keV band luminosity to bolomet- 
ric luminosity is weakly dependent on redshift and close 
to Unear in temperature (for kT > IkeV), means that a 
power law in mass and band luminosity implies a power 
law in mass and bolometric luminosity (and vice versa) 
to good approximation over the mass range of our data 
(> 3 X 10" Mq; see Figure 4). Performing this conversion, 
our best-fitting band luminosity-mass relation corresponds 
to a bolometric luminosity-mass relation with normalization 
0^"" = 1.26 ± 0.17 and slope /Jf™ = 1.59 ± 0.09. 

Finally, we note that our temperature-mass relation, 
determined with temperatures measured in the (0.15-l)r5oo 
annulus, is not significantly different from that obtained us- 
ing temperatures measured within the full region r < 7-500. 
For this center-included (ci) temperature-mass relation, we 
find normalization Z?^;^; = 0.90 ± 0.04, slope /S'^^i = 0.46 ± 
0.04, and marginal scatter (T(J„ = 0.059 ± 0.007 from the 
analysis of the XLF data. 



5.1.2 Results incorporating external cosmological data 

The precision of our results can be improved somewhat by 
incorporating the external cosmological data sets summa- 
rized in Section 2.4. These additional data reduce the pos- 
terior uncertainty on the cosmological parameters of the 
model, allowing the XLF data to more effectively constrain 
the scaling relations. (Recall that the parameters describing 
the scaling relations and the underlying cosmology are deter- 
mined simultaneously.) The results of this combined analysis 
are shown in Tables 7 and 8, and in Figure 2. Improvements 
are limited to the normalizations of the two relations, /3o'" 
and /3q^, and the luminosity-mass slope, /Ji"". The former 
are primarily degenerate with /gas, which sets the overall 



^■^ Wc refer specifically to the linmix_err routine, which does not 
account for selection bias. Since our own analysis is also Baycsian, 
it can be thought of as a specialization of the general approach 
advocated by Kelly (2007). 




Figure 4. Best-fitting luminosity-mass relations, where lumi- 
nosity is either bolometric or in the 0.1-2.4 keV band. The dot- 
dashed, blue line and dashed, red line show respectively the best- 
fitting bolometric and band luminosity relations, extrapolated to 
lower masses. The solid, black line shows the predicted band lu- 
minosity relation, assuming the best-fitting bolometric relation is 
valid at all masses and using the best-fitting temperature-mass 
relation. If the true, fundamental bolometric relation is a simple 
power law, a break in the slope of the band luminosity relation is 
expected at ~ 10^* Mq ; however, the band luminosity relation is 
still well described by a power law for A/ > 2 X lO^"* A/q . 



Table 8. Marginalized 68.3 per cent confidence intervals on the 
parameters of the temperature-luminosity relation from XLF 
data alone, and from the combination of all the cosmological 
data sets (Section 2.4). Tildes indicate parameters describing 
a bolometric (as opposed to 0.1-2. 4keV band) luminosity rela- 
tion, which are derived from the 0.1-2.4 keV results using the 
temperature-mass relation (Table 7). 



Parameter 



XLF only 



all data 



t-£ normalization 




0.55 


± 


0.08 


0.57 


± 


0.07 


t-£ slope 




0.30 


± 


0.08 


0.29 


± 


0.08 


t—£ scatter 




0.08 


± 


0.02 


0.08 


± 


0.02 


i-^bol normalization 


Pi' 


0.48 


± 


0.08 


0.50 


± 


0.07 


t-^bol slope 




0.27 


± 


0.04 


0.27 


± 


0.04 



mass scale, while the latter is degenerate with erg and influ- 
ences the total number of clusters that should be detected. 



5.1.3 The derived temperature-luminosity relation 

Given the model defined in Section 3, the nominal 
temperature-luminosity relation and its scatter are not triv- 
ially related to the fitted parameters, even at 2: = 0. The 
stochastic relation is 



P{l,t) _ J dmP{£,t\m)P{m) 



Pii) J dmP[l\m)P{m) 



(8) 



where P{l,t\m) is the bivariate scaling relation, P{l\m) is 
the marginal luminosity-mass relation, and P(m) is pro- 
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portional to the mass function. Thus, a detailed prediction 
for the temperature-luminosity relation requires specifying 
a set of cosmological parameters. Defining the nominal re- 
lation 



(9) 



with normal intrinsic scatter att, we find /3o ~ 0.55 ± 0.08, 
/3f = 0.30±0.08 and cr« = 0.08±0.02. Since the dependence 
of these estimates on the details of the mass function is not 
very strong, we have not determined these best fitting values 
and uncertainties precisely by marginalizing over cosmolog- 
ical parameters, although this would be straightforward in 
principle; instead, we have taken the less computationally 
intensive approach of using the mass function for our ref- 
erence cosmology at 2; = 0, and including in the error bars 
the difference between estimates based on the correspond- 
ing reference P(m) and a uniform P(m) as a conservative 
systematic allowance. These results, and the corresponding 
constraints for the bolometric temperature-luminosity rela- 
tion appear in Table 8. 



5.1.^ The derived Yx-mass relation 

The quantity Yx ~ M^askTce, an approximation to 
the thermal energy in the intracluster medium, has 
been proposed as a low scatter proxy for total mass 
(Kravtsov, Vikhlinin, & Nagai 2006), and has found appli- 
cation in cosmological work (Vikhlinin et al. 2009a, b). Using 
this definition, and the expectation that gas mass and to- 
tal mass are directly proportional (Mgas ~ /gasA/, with /gas 
constant), the virial theorem predicts the Yx-mass relation 
to be (cf. Equation 17) 



Yx (X M [E{z)Mf^^ 
Defining 

£(2:) Mgas (rsoo ) fcTcc 



y = logi 



101^ M<7, keV 



(10) 



(11) 



and allowing a general temperature-mass slope, /Si™, we 
can write the self-similarly evolving Fx-niass relation in 
terms of the gas mass fraction and the parameters of the 
temperature-mass relation: 

y = PE'^ + Pr^^ (12) 

= /So"" -I- logy, [/gas(r50o)l + {PT + 1) m. 

The best fit temperature-mass relation from the XLF data 
and mean /gas value at rsoo thus predict a normalization of 
/3^'" = 0.01 ± 0.05 and a slope of Z?^"" = 1.48 ± 0.04 for the 
Yx-mass relation (Table 9). 

This prediction is compared with the data in Figure 5. 
The scatter about the line is similar to the intrinsic scatter 
in the temperature-mass relation (~ 0.05, or 12 per cent), 
as expected from Equation 12 and the fact that intrinsic 
scatter in /gas for massive clusters remains undetected in 
current data (with a 68.3 per cent confidence upper bound 
< 7 per cent; A08). These considerations, and the fact that 



We note that the common practice of, e.g., estimating the 
temperature- luminosity slope as corresponds to the 

case of a fiat mass function, and is therefore biased. 



Table 9. Marginalized 68.3 per cent confidence intervals on the 
parameters of the Yx—raass relation, derived from the best-fitting 
temperature-mass relation and mean /gas value. The intrinsic 
scatter listed here is estimated from the scatter of the data about 
the derived relation, and is consistent with the value expected 
based on the scatter in the temperature— mass relation. 



Parameter 




Constraint 


y-m normalization 




0.01 ±0.05 


y-m slope 




1.48 ±0.04 


y-m scatter 




0.052 ±0.008 




(10^^ M^) 



Figure 5. Total mass and Yx = Afgas^Tcc are plotted for our 
best-fitting cosmology (Paper I) and compared with the Yx-mass 
relation predicted by the best-fitting temperature-mass relation 
and /gas value. The intrinsic scatter, indicated by the orange point 
with error bars, is similar to that of the temperature— mass rela- 
tion (~ 0.05). Note that the X and Y error bars are correlated, 
since Mgas is used to compute both Yx and the total mass, M500. 
To illustrate this correlation (p si 0.85), the (joint, 68.3 per cent 
confidence) error ellipse is shown in blue for the rightmost cluster. 



the gas mass can be measured both more precisely and with 
significantly smaller systematic uncertainties than tempera- 
ture^® for a given exposure time, suggest that Mgas may be 
a preferable mass proxy to either kT or Yx , at least for the 
most massive clusters where /gas is constant. 

5.2 Is the minimal model sufficient? 

An important question to address is whether the minimal 
model, as we have defined it, provides an acceptable de- 
scription of the data. Given the complexity and Bayesian 
nature of the model, testing the acceptability of the fit is 
necessarily more involved than following a simple Pre- 
scription. Our approach, which is not a unique solution to 
the problem, is to compare each of the following aspects of 
each cluster sample to the predictions of the best fit: 

^® Uncertainties in background modeling are currently the largest 
source of systematic error in temperature measurements at rsoo. 
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(i) the number of clusters detected, 

(ii) the two-dimensional distribution of clusters in red- 
shift and survey flux, 

(iii) the marginal distribution of luminosities from follow- 
up observations given the mass and redshift measurements, 

(iv) the marginal distribution of temperatures from 
follow-up observations given the measurements of mass, lu- 
minosity and redshift. 

The number of clusters in each sample is straightfor- 
wardly Poisson distributed. For each sample, we define a 
significance value as the probability of finding a number of 
clusters that is at least as discrepant from the predicted 
mean as the measured value. 

To test the survey flux-redshift distribution of sources, 
we use a two-dimensional analog of the Kolmogorov-Smirnov 
(K-S) test due to Peacock (1983) and Fasano & Franceschini 
(1987, see also Press et al. 1992). Briefly, the test involves 
comparing the predicted and detected numbers of clusters 
in each of the four quadrants {z ^ z, F ^ F; z < z, F ^ F; 
etc.) about each cluster detection, {z,F). The maximum 
discrepancy between these predictions and measurements is 
compared with a modification of the Kolmogorov D distri- 
bution to obtain the associated significance. 

The goodness of fit of the marginal luminosity-mass 
relation is tested by comparing the luminosities measured 
in the follow-up observations, £, to their predicted distri- 
bution given the measured masses and redshifts, m and z, 
and the detection of the cluster, indicated by /. This den- 
sity, P{i'\z,'m, I), is closely related to the probability as- 
sociated with a detected cluster in the likelihood function 
(see Appendix B), the principle difference being that the 
survey flux measurement should not be taken into account 
here. That is, the predicted distribution of £ accounts for 
the meiss function, the measured mass with its statistical 
error, and the selection bias. However, it should not ac- 
count for the specific flux value measured in the survey, 
or the expression would reduce to a test of how well the 
two flux measurements agree. This distribution can be cal- 
culated numerically as a function of £' and compared with 
the measured values, £. Since the shape of P{i'\z,m, I) is 
slightly different for each cluster, these residuals are not in- 
dividually statistically meaningful. However, the cumulative 

values C = P{i'\z,m,I)di' will be distributed uniformly 
on [0, 1], provided that the predicted distributions are ac- 
curate, independent of the fact that those distributions are 
not identical. Thus, we can test whether the fit is a consis- 
tent description of the luminosity-mass relation by checking 
whether these C values are uniformly distributed, using a 
one-sample K-S test. An analogous procedure is used to test 
the residuals of the temperature-mass relation. 

For the purpose of this goodness of fit test, we use the 
set of parameters producing the highest-likelihood sample 
in the Markov chains; this need not be precisely the mode 
of the posterior distribution, but should be adequate for the 
purpose of testing the fit. The 12 significance values (4 tests 
for each of the 3 cluster samples) obtained through the pro- 
cedure described here are displayed in Table 10. None of 
the individual tests produced a sufficiently low significance 
value (< 0.01) to indicate an inadequacy in the model. We 
note that if, for example, one test had produced a very small 
significance, it would still need to be interpreted in the con- 



Table 10. Significance values for each goodness of fit test applied 
to each cluster sample. Very low values (< 0.01) would indicate 
that the model is inadequate to describe the data (though, for- 
mally, the probability of randomly obtaining some low values in 
a sample of 12 would need to be taken into account). 



test 


BCS 


REFLEX 


MACS 


number 


0.73 


0.69 


0.21 


redshift-flux 


0.18 


0.78 


0.79 


luminosity-mass 


0.42 


0.93 


0.35 


temperature-mass 


0.19 


0.28 


0.92 



text of the larger goodness of fit test; that is, the likelihood 
of randomly obtaining one low significance value out of 12 
would need to be taken into account. 

5.3 Extensions to the model 

Although the minimal set of parameters produces an accept- 
able fit, we can still investigate whether any of the extensions 
to the simple model, summarized in Table 11, are preferred 
by the data. To do this, we have performed additional anal- 
yses with each of the extension parameters individually free. 
The additional data described in Section 2.4 were included in 
these analyses in order to constrain the cosmology as much 
as possible, maximizing our sensitivity to the scaling rela- 
tions. For the parameters controlling evolution with redshift, 
we choose the particular form of evolution Q{z) = 1 -I- z in 
Equations 6 and 7; however, our conclusions are identical if 
the other commonly chosen function, C,[z) = E{z), is used. 

The parameters controlling evolution in the marginal 
temperature-mass scatter and the scatter correlation coeffi- 
cient, a't^ and Pitmi were not constrained within the allowed 
region (priors in Table 11). The best-fitting temperature- 
mass scatter, atm = 0.055, is significantly smaller than the 
average measurement uncertainty on t, ~ 0.1, so the inabil- 
ity of the data to constrain evolution in the scatter is per- 
haps not surprising. Similarly, constraining the scatter cor- 
relation as a function of redshift would require many more 
data points, so that the correlation in £ and t at similar 
masses could be estimated at multiple redshifts. The rela- 
tively large uncertainty on pum (±0.19) indicates that this 
is a challenging measurement to make with the current data 
set, even without allowing evolution with redshift. 

The posterior distribution of the shape parameter of 
the luminosity-mass marginal scatter was sharply bounded 
at Xcrn. > —5, but was otherwise flat within the allowed 
region. Our inability to signiflcantly constrain asymmetry 
in the luminosity-mass scatter is consistent with the idea 
that the current data consist predominantly of clusters that 
are at least of average brightness given their masses. That is, 
the data do not sample well the lower tail of the distribution 
of luminosity given mass; hence, we can only rule out very 
negative skewness, where scatter to higher luminosities is 
suppressed beyond what is observed. 

The remaining parameters, governing the evolu- 
tion of the normalizations of the luminosity-mass and 
temperature-mass relations and of the luminosity-mass 
marginal scatter (Z?!™, /^I™ a-nd o-^„), are well constrained; 
the 68.3 per cent confidence intervals obtained for each are 
listed in Table 11. The fact that these constraints exclude 
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Table 11. Parameters which can extend the scahng relation model. Constraints of the form X±ax are 68.3 per cent confidence intervals; 
the bound on Af,„ refers to a sharp cut-off in posterior density. The notation U{a, b) indicates the uniform distribution with endpoints 
a and b. Deviance information criteria are with respect to the best fit using the minimal scaling relation model. 



Parameter 




Prior 


Constraint 


ADIC 


Extra evolution in £-m normalization 




W{-2,2) 


-0.5 ±0.4 


-0.9 


Extra evolution in t-m normalization 


P2 


W{-2,2) 


-0.2 ±0.2 


-1.6 


Evolution in £-m marginal scatter 




U{-2,2) 


-0.6 ±0.3 


-1.2 


Evolution in t-m marginal scatter 




U(-2,2) 


none 




Evolution in £-t scatter correlation 


Petm 


U(-2,2) 


none 




Asymmetry in £—m marginal scatter 




ul-W, 10) 


> -5 





zero at the l-2a level is not itself sufficient to conclude that 
the data require additional parameters in the model. In- 
stead, a model selection test must be employed to determine 
whether the observed improvement in likelihood is signifi- 
cant compared to that expected to occur randomly when a 
new parameter is added. We address this question using the 
Deviance Information Criterion (DIG) of Spiegelhalter et al. 
(2002, see also Liddle 2007). This information criterion is 
well suited to Bayesian problems and, unlike the Bayesian 
evidence, does not depend strongly on the somewhat arbi- 
trary width of the priors. The DIG is 

DIG = 21n£ - 2inZ (13) 

where £, is the likelihood, the bar indicates an average over 
the posterior, and £. is the likelihood evaluated at some mea- 
sure of center of the posterior such as the mean or mode. We 
use the maximum likelihood sampled by the Markov chains 
as £. 

The difference in DIG obtained for these three tests 
compared with that from the minimal model is listed in Ta- 
ble 11. The ADIG values range from —0.9 to —1.6, which, 
according to the Jeffreys' scale conventionally used to inter- 
pret these values, indicates that none of the extensions to 
the model is strongly favored. We conclude that the mini- 
mal scaling relation model provides an adequate description 
of the data, and conversely that the data do not prefer more 
complex models at a statistically significant level. 

5.4 Comparison with other work 

Previous results on the slopes and scatters of the scaling rela- 
tions studied here have varied widely. Since no other authors 
have accounted fully for the underlying cluster mass function 
and the sample selection function, as we do here, such large 
systematic variations could be expected to result from the 
use of different data sets and fitting methods. Moreover, our 
analysis uses Mgas as a proxy for total mass, which should 
be more reliable in the mass range studied here than the hy- 
drostatic mass estimates used in most previous work. This 
is particularly true for the temperature-mass relation, since 
hydrostatic mass measurements are highly correlated with 
measured temperatures. In addition, the 21 January 2009 
Chandra calibration update may be significant; the reduc- 
tion of relative effective area at high energies disproportion- 
ally reduces temperatures of hot clusters compared with cool 

^'^ Typically, a ADIC of < —5 is interpreted as "strong" evidence 
for the more complex model, and < —10 is "conclusive." 



systems, resulting in a net fiattening of the temperature- 
mass relation. We note also the possibility that the slopes 
of the scaling relations are a function of mass, in which case 
our results should only be compared with those of similarly 
high-mass cluster samples. Nevertheless, for completeness, 
we compare in this section our constraints with other recent 
results in the literature. 

For the luminosity-mass relation, we find a slope of 
/3f™ = 1.29±0.07 for 0.1-2.4 keV luminosities, or 1.59±0.09 
for bolometric luminosities. This slope is steeper than that 
expected from the virial theorem (4/3 for the bolometric 
slope), but shallower than most other results, which typi- 
cally span the range 1.5-1.8 for band luminosities or 1.7-2.3 
for bolometric luminosities (e.g. Reiprich & Bohringer 2002; 
Zhang et al. 2007; M08; Rykoff et al. 2008; Zhang et al. 
2008; Pratt et al. 2009; Vikhfinin et al. 2009a), though we 
note good agreement with a previous study of very lumi- 
nous clusters (Allen et al. 2003). 

Our temperature-mass slope, Pl"^ = 0.48 ± 0.04, is 
shallower than the virial prediction of 2/3 and most oth- 
ers in the literature (typically 0.6-0.7, e.g. Allen et al. 2001; 
Finoguenov et al. 2001; Arnaud et al. 2005; Popesso et al. 
2005; Vikhlinin et al. 2006; Sun et al. 2009), though there 
do exist some results with similar slopes (e.g. Ettori et al. 
2004). As noted above, some of the differences may be due 
to our use of gas mass as a proxy for total mass, which 
should be more reliable than hydrostatic masses for studying 
the temperature-mass relation. Consistently, our estimate 
of the Yx-mass slope is also somewhat shallower than that 
derived in other studies (e.g. Arnaud et al. 2007; Sun et al. 
2009; Vikhlinin et al. 2009a). 

Some authors studying a combination of low and high 
redshift clusters have remarked that the evolution expected 
from self similarity appears empirically justified, in that low 
and high redshift clusters with the appropriate E{z) factors 
appear to lie about the same line (e.g. Maughan et al. 2008; 
Zhang et al. 2008). Ettori et al. (2004) and Morandi et al. 
(2007) went further, by adopting a model for evolution as in 
Equation 6 with = 1+z. Using a modified least-squares fit, 
they found no evidence for evolution beyond the self-similar 
expectation, in agreement with our conclusions. 

Pratt et al. (2009) and Vikhlinin et al. (2009a) found 
luminosity-mass (0. 1-2.4 keV band) slopes ~ 1.6 using a 
procedure in which the luminosity data are corrected for 
the expected Malmquist bias, given a flux limit or full selec- 
tion function. Vikhlinin et al. (2009a) also considered the 
possibility of departures from self-similar evolution, find- 
ing a marginal preference for non-self-similar evolution of 
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E(z)M,„„(10^^MJ 

Figure 6. Luminosity-mass data for our reference eosmology, 
where luminosities include only emission from 0.15 < r/rsoo < 
1. Note that the X and Y error bars are anti-correlated (p « 
—0.4), as illustrated by the blue (joint, 68.3 per cent confidence) 
error ellipse shown for the leftmost cluster; this anti-correlation 
is accounted for in the analysis. Black points indicate data used 
in our standard analysis; blue points are high-redshift (0.5 < z < 
0.7) Faint MACS clusters from Ebeling et al. (2007), and green 
points are 400d clusters (0.35 < z < 0.9; Burenin et al. 2007). 
The best fit to this center-excised relation, shown in red, has a 
slope in agreement with the virial theorem and is consistent with 
having intrinsic scatter at the few per cent level. The data are 
consistent with self-similar evolution. 

the luminosity-mass relation. Their estimate of the intrinsic 
scatter in the relation agrees well with our own. 



6 CENTER-EXCISED LUMINOSITY SCALING 
RELATIONS 

It has long been known that the luminosity of clusters har- 
boring cooling cores is significantly enhanced by the emis- 
sion from bright, cool gas in their centers (e.g. Fabian et al. 
1994; Allen & Fabian 1998; Markevitch 1998; Peres et al. 
1998). More recently, a dramatic reduction in luminosity- 
mass scatter has been demonstrated when luminosities are 
measured excluding cluster centers (typically r < O.lSrsoo; 
Maughan 2007; Zhang et al. 2007). Our data confirm and 
reinforce the conclusion that the center-excised luminosity- 
mass relation is significantly tighter than the center-included 
relation, as Figure 6 shows. (Note that, for clarity, error bars 
in the figure do not include systematic uncertainties that are 
used elsewhere in our analysis; see below.) 

To the extent that the center-excised relation exhibits 
zero or negligible intrinsic scatter, selection effects will not 
bias the results of regression performed on these data. [How- 
ever, it is still important to marginalize over cosmological 
parameters when fitting the data, since they are reflected 
in the values of E{z).] If selection bias is not an issue, the 
strict requirements used to deflne the data set can be re- 
laxed; Figure 6 shows the standard data set used through- 
out this work in black, with an additional 12 high-redshift 



(0.5 < 2 < 0.7) Faint MACS clusters from Ebeling et al. 
(2007) in blue, and 22 clusters from the 400 Square De- 
gree survey (400d, 0.35 < z < 0.9; Burenin et al. 2007) in 
green. The high-redshift MACS and 400d data satisfy Equa- 
tion 2 and the respective flux limits 10~^^ ergs~^ cm~^ (0.1- 
2.4 keV band) and 1.4 x 10"" ergs"^ cm"^ (0.5-2.0 keV)." 
Chandra data for these clusters were reduced as described 
in Section 2; details and results of the observations appear 
in Tables 13 and 14. 
We flt the relation 

(4c (m)) = 7?" + + 7?" log.o (C) , (14) 

where ica is deflned as in Equation 5, but with luminos- 
ity from within radius 0.15r5oo of the center excluded (see 
Section 2.2.3). The fit was performed using the Bayesian 
linear mixture method of Kelly (2007), which accounts for 
the correlation in mass and luminosity measurement er- 
rors.^® Somewhat counterintuitively, the correlation between 
center-excised luminosity and mass measurements is neg- 
ative; although both quantities depend on the measured 
electron density, a stronger anti-correlation arises because 
a larger mass implies a larger rsoo, which results in a net 
reduction in luminosity due to the center-excision. 

The fit was marginalized over the cosmological param- 
eter values allowed by the combination of XLF and other 
cosmological data (Paper I). Note that, unlike elsewhere in 
this paper, we do not simultaneously fit for cosmological 
parameters here, but only marginalize the results over the 
allowed set of cosmological models; hence the fits in this sec- 
tion directly use only the clusters shown in Figure 6, not the 
additional clusters that were detected in the RASS but not 
followed up. 

From the standard data set (BCS, REFLEX and Bright 
MACS), we find normalization 70"" = 0.71 ± 0.06 and slope 
7f'" = 1.02 ± 0.05, wfth an upper bound of 0.015 on the 
intrinsic scatter (68.3 per cent confidence), for the self- 
similarly evolving (72™ = 0) case. The extended data set, 
using the high-redshift MACS clusters and 400d clusters, 
yields virtually identical constraints. Similar results are ob- 
tained with the other simple fitting procedures discussed in 
Section 5.1. These results, and those for the corresponding 
bolometric luminosity-mass relation, are shown in Table 12. 
Note that, unlike the procedure in Section 5.1, fitting the 
bolometric relation in this case does not require the use of 
a temperature-mass relation; we directly fit the masses and 
bolometric luminosities derived from the follow-up X-ray ob- 
servations. 

We stress that the limit on the intrinsic scatter above 



For the 400d clusters, we also require the core-included 0.1- 
2.4 keV luminosity from follow-up observations to be > 2 X 
10** erg s~^ , as the few clusters that otherwise meet our selection 
criteria have extremely poor signal-to-noise. The high-redshift 
MACS clusters are all above this luminosity limit. 

We note that it is extremely straightforward, within the 
Bayesian framework, to account for correlation in the mass and 
luminosity measurements. In brief, the likelihood function inte- 
grates over the true values of mass and luminosity for every clus- 
ter, weighted by the sampling distribution, which gives the proba- 
bility of obtaining the measured values given the true values. The 
sampling distribution in this calculation can easily be bivariate, 
with non-zero correlation. See also Gelman et al. (2004). 
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Table 13. Details of the follow-up observations of high-redshift (2 > 0.5) MACS clusters and 400d clusters (see caption for Table 1). 



Name RA (J2000) Dec (J2000) Date Detector Mode Exposure (ks) 



MACS J2214.9-1359 


22 14 57.3 


-14 


00 


11 


2002 


Dec 


22 


ACIS- 


-I 


VFAINT 


14.6 
















2003 


Nov 


17 


ACIS- 


-I 


VFAINT 


14.7 


MACS J0911. 2+1746 


09 11 


10.9 


+17 


46 


31 


2003 


Feb 


23 


ACIS- 


-I 


VFAINT 


17.4 
















2004 


Mar 


08 


ACIS- 


-I 


VFAINT 


23.8 


MACS J0257.1-2325 


02 57 9.1 


-23 


26 


04 


2000 


Oct 


03 


ACIS- 


-I 


FAINT 


19.6 
















2003 


Aug 


23 


ACIS- 


-I 


VFAINT 


18.2 


MACS J0454.1-0300 




04 


114 
±1.4 


-03 


00 


51 


2000 


Jan 


14 


ACIS- 


-I 


VFAINT 


13.8 
















2000 


Oct 


08 


ACIS- 


-s 


FAINT 


39.9 


MACS J1423.8+2404 


1 A 


Zo 


/IV 

4 ( .y 


+24 


04 


43 


2001 


Jun 


01 


ACIS- 


-I 


VFAINT 


18.5 
















2003 


Aug 


18 


ACIS- 


-s 


VFAINT 


114.5 


MACS J1149.5+2223 


1 1 
± 1 


AQ 


0O.4 


+22 


24 


04 


2001 


Jun 


01 


ACIS- 


-I 


VFAINT 


17.0 
















2003 


Feb 


07 


ACIS- 


-I 


VFAINT 


17.7 


MACS J0717.5+3745 


(17 


1 7 


32 1 


+37 


45 


21 


2001 


Jan 


29 


ACIS- 


-I 


FAINT 


16.5 
















2003 


Jan 


08 


ACIS- 


-I 


VFAINT 


55.1 


MACS J0018.5+1626 


on 


1 S 


00.4 


+16 


26 


13 


2000 


Aug 


18 


ACIS- 


-I 


VFAINT 


66.9 


MACS J0025.4-1222 


nn 
uu 


ZO 


9Q Q 

zy .y 


-12 


22 


45 


2002 


Nov 


11 


ACIS- 


-I 


VFAINT 


15.7 
















2004 


Aug 


09 


ACIS- 


-I 


VFAINT 


24.4 


MACS J2129.4-0741 


91 


9Q 


9^^ 7 
ZO. / 


-07 


41 


31 


2002 


Dec 


23 


ACIS- 


-I 


VFAINT 


10.5 
















2003 


Oct 


18 


ACIS- 


-I 


VFAINT 


19.3 


MACS J0647.7+7015 


UD 


A7 
4 1 


AQ 7 
4y. i 


+70 


14 


56 


2002 


Oct 


31 


ACIS- 


-I 


VFAINT 


16.7 
















2003 


Oct 


07 


ACIS- 


-I 


VFAINT 


19.9 


MACS J0744.8+3927 


u / 


AA 


•^9 
DZ.o 


+39 


27 


27 


2001 


Nov 


12 


ACIS- 


-I 


VFAINT 


19.7 
















2003 


Jan 


04 


ACIS- 


-I 


VFAINT 


15.9 
















2004 


Dec 


03 


ACIS- 


-I 


VFAINT 


49.5 


400d J0302. 3-0423 


uo 


(19 
uz 


21 1 


-04 


23 


24 


2005 


Dec 


07 


ACIS- 


-I 


VFAINT 


9.0 


400d J1212. 3+2733 


12 


12 


18 4 


+27 


33 


02 


2005 


Mar 


17 


ACIS- 


-I 


VFAINT 


12.5 


400d J0318. 5-0302 


uo 


1 s 


00.0 


-03 


02 


58 


2005 


Mar 


15 


ACIS- 


-I 


VFAINT 


10.7 


400d J0809. 6+2811 


uo 


(1Q 

uy 


41 (1 
41 .U 


+28 


12 


01 


2004 


Nov 


30 


ACIS- 


-I 


VFAINT 


5.3 


400d J1416. 4+4446 




ID 


9fi 1 
Zo. 1 


+44 


46 


43 


1999 


Dec 


02 


ACIS- 


-I 


VFAINT 


25.0 


400d J1701. 3+6414 


1 7 


(11 
Ul 


9*^ 1 
Zo. 1 


+64 


14 


08 


2000 


Oct 


31 


ACIS- 


-I 


VFAINT 


37.1 


400d J0355. 9-3741 




00 


oy .4 


-37 


41 


46 


2006 


Jan 


12 


ACIS- 


-I 


VFAINT 


22.3 


400d J1002. 1+6858 


1 n 
±u 


(19 

uz 


(1Q (1 

uy .u 


+68 


58 


36 


2005 


Jan 


05 


ACIS- 


-I 


VFAINT 


17.9 


400d J0030.5+2618 


no 
uu 


ou 


00.0 


+26 


18 


09 


2005 


May 


28 


ACIS- 


-I 


VFAINT 


11.8 


400d J1524.6+0957 


1 

±0 


9A 
Z4 


oy .0 


+09 


57 


46 


2002 


Apr 


01 


ACIS- 


-I 


VFAINT 


40.7 


400d J1357.3+6232 


1 


1 


1 7 
1 1 .u 


+62 


32 


51 


2006 


Jan 


24 


ACIS- 


-I 


VFAINT 


21.0 
















2006 


Jan 


29 


ACIS- 


-I 


VFAINT 


15.8 


400d J1120.9+2326 


11 


20 


57.5 


+23 


26 


29 


2001 


Apr 


23 


ACIS- 


-I 


VFAINT 


55.0 


400d J0956.0+4107 


09 


56 


03.3 


+41 


07 


08 


2003 


Dec 


30 


ACIS- 


-I 


VFAINT 


15.0 
















2005 


Jan 


28 


ACIS- 


-I 


VFAINT 


35.6 


400d J0328. 6-2140 


03 


28 


36.2 


-21 


40 


22 


2005 


Mar 


15 


ACIS- 


-I 


VFAINT 


29.4 
















2005 


Mar 


18 


ACIS- 


-I 


VFAINT 


3.7 


400d J1120.1+4318 


11 


20 


06.9 


+43 


18 


05 


2005 


Jan 


11 


ACIS- 


-I 


VFAINT 


15.2 


400d J1334.3+5031 


13 


34 


20.1 


+50 


31 


01 


2005 


Aug 


05 


ACIS- 


-I 


VFAINT 


7.0 


400d J0542.8-4100 


05 


42 


50.1 


-41 


00 


04 


2000 


Jul 


26 


ACIS- 


-I 


FAINT 


41.1 


400d J1202. 2+5751 


12 


02 


17.8 


+57 


51 


54 


2005 


Sep 


02 


ACIS- 


-I 


VFAINT 


45.5 


400d J0405.4-4100 


04 


05 


24.5 


-41 


00 


19 


2005 


Oct 


27 


ACIS- 


-I 


VFAINT 


6.6 
















2006 


May 


19 


ACIS- 


-I 


VFAINT 


57.7 


400d J1221. 3+4918 


12 


21 


26.1 


+49 


18 


31 


2001 


Aug 


05 


ACIS- 


-I 


VFAINT 


62.7 


400d J0152. 7-1358 


01 


52 


41.1 


-13 


58 


07 


2000 


Sep 


08 


ACIS- 


-I 


FAINT 


29.2 


400d J1226. 9+3332 


12 


26 


57.9 


+33 


32 


49 


2000 


Jul 


31 


ACIS- 


-s 


VFAINT 


9.8 
















2003 


Jan 


27 


ACIS- 


-I 


VFAINT 


26.1 
















2004 


Aug 


07 


ACIS- 


-I 


VFAINT 


25.3 



must be treated with caution. The luminosities and masses 
determined from follow-up observations incorporate various 
systematic allowances, most importantly those accounting 
for instrument calibration and the cluster-to-cluster varia- 
tion of /gas (Section 2.2.3). The precise definition of these 
allowances has very little impact on the mean relation de- 
rived here, or on the results of Section 5; however, because of 
its small size, the scatter in the center-excised luminosity- 



mass relation is sensitive to these allowances. To estimate 
the effect that the systematic allowances have on the de- 
termination of the intrinsic scatter, we repeat the fit using 
statistical error bars only (as shown in Figure 6). In this 
case, we find an intrinsic scatter of 0.025 ± 0.004, or ~ 6 
per cent. We conclude that, while the definition of a formal 
confidence interval for the center-excised intrinsic scatter is 
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Table 14. Redshifts and derived properties of high-rcdshift (z > 0.5) MACS clusters and 400d clusters from follow-up observations (see 
caption for Table 4). Note that these data are used only in Section 6. 



Name z rsoo A/gas,500 M500 L^o^cc 

(Mpc) (10^'* A/©) (10" M©) (lO^^ergs-i) 



MACS J2214.9-1359 


0.5027 


1.39 


± 


0.08 


1.51 


± 


0, 


.22 


13.2 


± 


2.3 


8.28 


± 


0.37 


MACS J0911. 2+1746 


0.5049 


1.22 


± 


0.06 


1.03 


± 


0, 


.11 


9.0 


± 


1.2 


5.41 


± 


0.25 


MACS J0257.1-2325 


0.5049 


1.20 


± 


0.06 


0.97 


± 


0, 


.13 


8.5 


± 


1.3 


6.59 


± 


0.30 


MACS J0454.1-0300 


0.5377 


1.31 


± 


0.06 


1.31 


± 


0, 


.13 


11.5 


± 


1.5 


9.50 


± 


0.41 


MACS J1423.8+2404 


0.5431 


1.09 




0.05 


0.76 




0, 


.08 


6.6 


± 


0.9 


4.64 


± 


0.20 


MACS J1149. 5+2223 


0.5444 


1.53 


± 


0.08 


2.13 


± 


0, 


.28 


18.7 


± 


3.0 


12.61 




0.52 


MACS J0717.5+3745 


0.5458 


1.69 


± 


0.06 


2.85 


± 


0, 


,24 


24.9 


± 


2.7 


18.03 


± 


0.71 


MACS J0018.5+1626 


0.5456 


1.47 


± 


0.08 


1.88 


± 


0, 


.24 


16.5 


± 


2.5 


11.84 


± 


0.67 


MACS J0025.4-1222 


0.5843 


1.12 


± 


0.04 


0.87 


± 


0, 


.08 


7.6 


± 


0.9 


6.81 


± 


0.31 


MACS J2129. 4-0741 


0.5889 


1.25 


± 


0.06 


1.22 


± 


0, 


,14 


10.6 


± 


1.4 


7.88 


± 


0.39 


MACS J0647.7+7015 


0.5907 


1.26 


± 


0.06 


1.25 


± 


0, 


,15 


10.9 


± 


1.6 


8.07 


± 


0.37 


MACS J0744.8+3927 


0.6976 


1.26 


± 


0.06 


1.43 


± 


0, 


,15 


12.5 


± 


1.6 


9.07 


± 


0.34 


400d J0302.3-0423 


0.350 


0.99 


± 


0.07 


0.46 


± 


0, 


,09 


4.0 


± 


0.8 


2.52 


± 


0.24 


400d J1212.3+2733 


0.353 


1.35 


± 


0.10 


1.18 


± 


0, 


,22 


10.3 


± 


2.1 


5.20 


± 


0.35 


400d J0318. 5-0302 


0.370 


0.85 


± 


0.11 


0.30 


± 


0, 


,07 


2.6 


± 


0.7 


1.96 


± 


0.17 


400d J0809.6+2811 


0.399 


1.07 


± 


0.11 


0.62 


± 


0, 


,15 


5.4 


± 


1.4 


3.23 


± 


0.30 


400d J1416.4+4446 


0.400 


0.83 


± 


0.05 


0.29 


± 


0, 


,05 


2.5 


± 


0.5 


1.71 


± 


0.12 


400d J1701.3+6414 


0.453 


0.90 


± 


0.06 


0.38 


± 


0, 


,07 


3.4 


± 


0.7 


2.31 


± 


0.14 


400d J0355. 9-3741 


0.473 


0.83 


± 


0.05 


0.32 


± 


0, 


,06 


2.8 


± 


0.5 


1.65 


± 


0.15 


400d J0030.5+2618 


0.500 


0.82 


± 


0.05 


0.31 


± 


0, 


,06 


2.7 


± 


0.5 


1.96 


± 


0.21 


400d J1002. 1+6858 


0.500 


0.93 


± 


0.07 


0.46 


± 


0, 


,09 


4.0 


± 


0.8 


2.16 


± 


0.22 


400d J1524.6+0957 


0.516 


0.94 


± 


0.07 


0.49 


± 


0, 


,08 


4.2 


± 


0.8 


2.81 


± 


0.20 


400d J1357.3+6232 


0.525 


0.81 


± 


0.05 


0.31 


± 


0, 


,05 


2.7 


± 


0.5 


1.93 


± 


0.12 


400d J1120.9+2326 


0.562 


0.85 


± 


0.16 


0.39 


± 


0, 


,10 


3.4 


± 


0.9 


2.27 


± 


0.16 


400d J0956.0+4107 


0.587 


0.78 


± 


0.06 


0.30 


± 


0, 


,05 


2.6 


± 


0.4 


2.15 


± 


0.12 


400d J0328.6-2140 


0.590 


0.86 


± 


0.07 


0.40 


± 


0, 


,10 


3.5 


± 


0.9 


2.44 


± 


0.19 


400d J1120.1+4318 


0.600 


0.98 


± 


0.07 


0.61 


± 


0, 


,12 


5.3 


± 


1.1 


4.61 


± 


0.41 


400d J1334.3+5031 


0.620 


0.85 


± 


0.10 


0.41 


± 


0, 


,13 


3.6 


± 


1.2 


2.27 


± 


0.63 


400d J0542. 8-4100 


0.642 


0.99 


± 


0.08 


0.64 


± 


0, 


,14 


5.6 


± 


1.3 


3.64 


± 


0.26 


400d J1202. 2+5751 


0.677 


0.80 


± 


0.05 


0.36 


± 


0, 


,07 


3.2 


± 


0.7 


2.56 


± 


0.23 


400d J0405. 4-4100 


0.686 


0.72 


± 


0.04 


0.26 


± 


0, 


,04 


2.3 


± 


0.4 


2.36 


± 


0.16 


400d J1221.3+4918 


0.700 


0.97 


± 


0.05 


0.65 


± 


0, 


,09 


5.7 


± 


0.9 


4.27 


± 


0.27 


400d J0152. 7-1358 


0.833 


0.97 


± 


0.26 


0.89 


± 


0, 


,33 


7.8 


± 


3.0 


6.66 


± 


0.51 


400d J1226. 9+3332 


0.888 


1.00 


± 


0.05 


0.89 


± 


0, 


,10 


7.8 


± 


1.1 


7.66 


± 


0.33 



problematic, the scatter is likely to be very small (< 10 per 
cent). 

The small size of this scatter in the self-similar case 
strongly suggests that departures from self-similar evolution 
in the center-excised relation are also small. We test this pos- 
sibility by fitting Equation 14 with the additional parameter 
71" free, finding 7!™ = -0.12l° j^ from the standard data 
set and 72™ = —0.22 ± 0.25 from the extended data, both 
consistent with zero. The latter corresponds to a 10 per cent 
constraint on excess evolution in the normalization of the re- 
lation within z — 0.5 (or 16 per cent within 2 — 0.9). The 
constraints on the slope and scatter are similar to those in 
Table 12. 

The existence of such a simply evolving, low-scatter 
proxy for total mass is exciting in the context of future X- 
ray cluster searches. Since the center-excised luminosity at 
soft energies is derived from the count rate with little tem- 
perature dependence, it can be measured precisely without 
the need for high-quality spectral information. This is in 
contrast to other low-scatter mass proxies such as the aver- 
age temperature {kT) or thermal energy {Yx = MgasfcTcc), 
which are generally limited by the precision of temper- 



ature measurements. Although the upcoming eROSITA^" 
survey is unlikely to have sufficient spatial resolution to 
directly measure center-excised fluxes, a mission along the 
lines of the proposed Wide-Field X-ray Telescope (WFXT, 
Murray et al. 2008) could potentially take advantage of this 
tight relation to construct an effectively mass-limited clus- 
ter sample directly from an X-ray center-excised flux-limited 
survey. 

If the scatter between mass and center-excised lu- 
minosity is negligible, the dependence of the center- 
excised temperature-luminosity relation on the mass func- 
tion (Equation 8) vanishes. The parameters of the nominal, 
center-excised temperature-luminosity relation, 

(t(4c)> =7" + 7f4e, (15) 
can then be estimated simply as (Table 12) 

7o = /3o -Pi (16) 

ptm 
U _ Pi 

http : //www .mpe .mpg. de/projects .html#erosita 
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Table 12. Marginalized 68.3 per cent confidence intervals on 
power law models describing the center-excised luminosity-mass 
and temperature— luminosity relations from our standard and ex- 
tended data sets (see text). 70 and 71 are the normalization and 
slope of each relation when luminosities are measured in the 0.1- 
2.4 keV band; 70 and 71 describe the fits using bolomctric lu- 
minosities. Because the center excision effectively eliminates the 
intrinsic scatter in the luminosity-mass relation (residual system- 
atic scatter is bounded at the few per cent level in both cases), the 
data can be fit directly, without using the full apparatus discussed 
in Section 4 and Paper I. These constraints are marginalized over 
the range of cosmological parameter values allowed by the com- 
bination of XLF and other cosmological data. 



Parameter Standard Extended 



Icc-m normalization 


7r 


0.71 


± 


0.06 


0.71 


± 


0.06 


ice—'m slope 




1.02 


± 


0.05 


1.03 


± 


0.04 


t—ice normalization 




0.54 


± 


0.05 


0.55 


± 


0.05 


t-ice slope 




0.48 


± 


0.04 


0.47 


± 


0.04 


^bol cc^"^ normalization 


7r 


1.16 




0.08 


1.16 


± 


0.08 


4ol,cc-"i slope 




1.31 


± 


0.05 


1.30 


± 


0.05 


t~ibol ce normalization 




0.45 


± 


0.06 


0.45 


± 


0.06 


*-^bol,ce slope 


7f 


0.37 


± 


0.03 


0.38 


± 


0.03 



The intrinsic scatter is identical to that of the temperature- 
mass relation. 



cool cores (e.g. Voit & Bryan 2001; Peterson & Fabian 2006; 
McNamara & Nulsen 2007). Importantly, the consistency of 
our data with a simple, self-similarly evolving power-law 
model indicates that the effects of the excess heat on the 
global luminosity and temperature have not evolved strongly 
since redshift 0.5. 

Interestingly, when the central O.lSrsoo is excluded from 
the luminosity measurements, the slope of the bolometric 
luminosity-mass relation, 7^"" = 1.30 ± 0.05, is close to the 
virial prediction. This observation suggests that the excess 
heating is limited to the central regions of clusters, arguing 
in favor of local heating mechanisms, such as AGN feed- 
back, and against mechanisms that would heat the ICM 
more uniformly, such as preheating. The fact that the center- 
excised temperature-mass slope does not agree with the 
virial prediction does not necessarily change this picture, 
since that prediction applies to the mass-weighted temper- 
ature rather than the emission-weighted temperatures re- 
ported here.^'^ Note also that the emission-weighted tem- 
perature is more influenced by gas at small radii, i.e. closer 
to the source of excess heating, where the signal-to-noise is 
much higher than in the outskirts. The small intrinsic scat- 
ter of the center-excised luminosity-mass (particularly) and 
temperature-mass relations over a wide range in redshift, 
including terms for self-similar evolution as we have, fur- 
ther reinforces the notion that the self-similar model closely 
describes the ICM outside of cluster centers. 



7 ASTROPHYSICAL IMPLICATIONS 

The virial theorem makes specific predictions for the slopes 
of the scaling relations linking cluster mass, bolometric lu- 
minosity and mass-weighted temperature; namely 

Lhoi 



E{z) 



[E{z)M] 



4/3 



2/3 



[E{z)M] 

r 11/2 



(17) 



E(z) 



for quantities within a fixed critical-overdensity radius. 
Compared to this prediction, our result for the (center- 
included) bolometric luminosity-mass slope, /3f'" = 1.63 ± 
0.06, is significantly steeper, and our constraint on the 
temperature-mass slope, /Jj™ = 0.49 ± 0.04, is significantly 
shallower. These results are consistent with a picture in 
which the intracluster medium has been subject to excess 
heating;'^! energy injection heats the gas, raising the tem- 
perature and decreasing the density (hence the luminosity), 
with the effect being more prominent in lower mass systems. 
Condensation of the coldest gas in cluster cores, leaving be- 
hind hotter, less dense gas, may also play a role, although 
some additional heating is still required in order to account 
for the typically modest star formation rates observed in 



■^1 However, as noted below, the fact that emission-weighted 
rather than mass-weighted temperatures are measured compli- 
cates the physical interpretation of the temperature-mass slope. 
In this work, the luminosity— mass relation is preferred for study- 
ing the effects of excess heating of the ICM. 



8 CONCLUSION 

We have presented constraints on the X-ray scaling relations 
of galaxy clusters using a new analysis method, detailed 
in Paper I, that self-consistently accounts for all observa- 
tional biases and degeneracies with cosmological parame- 
ters. The data are drawn from three wide area. X-ray flux- 
limited cluster samples based on the RASS, span the red- 
shift range < z < 0.5, and satisfy the luminosity threshold 
L > 2.5 X 10**/i^Q^ erg s~^. Masses, X-ray luminosities, and 
average ICM temperatures at rsoo were derived from follow- 
up X-ray observations of 94 of the 238 flux-selected clus- 
ters, and used to constrain joint luminosity-temperature- 
mass scaling relations, marginalizing over a ACDM cosmo- 
logical background as well as conservative systematic al- 
lowances. Our analysis for the first time incorporates the 
cluster mass function and sample selection function in a fully 
self-consistent and statistically rigorous way. 

The data are consistent with a relatively simple picture 
in which the nominal temperature-mass and luminosity- 
mass relations in the observed mass range are described 
by power laws that evolve self similarly. The scatter in lu- 
minosity and temperature at fixed mass is adequately de- 
scribed by a simple, bivariate log-normal distribution. De- 
partures from self-similar evolution, evolution in the scatter, 
and asymmetry in the scatter were considered, but the data 
do not require any of these additions to the model at a sta- 
tistically significant level. Incorporating additional cosmo- 
logical data to minimize the degeneracies with cosmological 

Unfortunately, it is not currently feasible to measure mass- 
weighted temperatures within rsoo due to systematic uncertain- 
ties in the background modeling. 
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parameters, we find that the slope of the (center-included) 
luminosity-mass relation is = 1.34 ± 0.05 (correspond- 
ing to bolometric luminosity-mass slope 1.63 ± 0.06) while 
the slope of the (center-excised) temperature-mass relation 
is /Jj™ = 0.49 ± 0.04. The marginal scatters of the relations 
are atm = 0.185 ± 0.019 and at™ = 0.055 ± 0.008 (68.3 
per cent confidence intervals). The correlation coefficient 
of the luminosity and temperature scatters at fixed mass 
is relatively poorly constrained, but consistent with zero, 
pttm = 0.09±0.19. The resulting temperature-luminosity re- 
lation has slope /?" = 0.29 ±0.08 and intrinsic scatter an = 
0.08 ± 0.02. The steep slope of the bolometric luminosity- 
mass relation and shallow slope of the temperature-mass 
relation relative to the predictions of the virial theorem 
(respectively 4/3 and 2/3) are consistent with the now- 
standard picture of excess heating of the ICM. Moreover, 
the agreement of our data with self-similar evolution indi- 
cates that the heating mechanism was in operation before 
z — 0.5, and that its effects on global cluster properties do 
not evolve strongly at redshifts z < 0.5. 

Previous observations that the scatter in the 
luminosity-mass relation is dominated by the emission 
in the cluster center (r < 0.15r5oo) are reinforced by our 
data. Here, we find three main, new results. First, the 
intrinsic scatter in the center-excised luminosity-mass 
relation is consistent with zero; we derive an upper bound 
at the < 10 per cent level. Second, the center- excised, 
bolometric luminosity-mass slope, "f-f™ = 1.30 ± 0.05, is in 
agreement with the virial prediction. Third, we find that 
the evolution of this relation to be consistent with the 
self-similar expectation, and bound departures from self 
similarity at the 16 per cent level over the range z < 0.9 (10 
per cent over z < 0.5). Together, these results argue that 
the excess heating in clusters has been primarily limited to 
their centers, supporting local heating mechanisms such as 
AGN feedback over more global preheating. 

The existence of a tight, simply evolving correlation be- 
tween center-excised luminosity and mass is exciting for cos- 
mological applications, since the center-excised luminosity 
may provide a simpler and cheaper (in terms of required ex- 
posure time) mass proxy than average temperature or ther- 
mal energy. Potentially, it could form the basis for the se- 
lection of an effectively mass-limited cluster sample from a 
future X-ray survey, given an instrument with sufficient spa- 
tial resolution. 

In this paper, we have not presented any cosmological 
results, despite the fact that the analysis used here simulta- 
neously constrains cosmological parameters. Those results, 
as well as constraints on more complex cosmological mod- 
els, can be found in Paper I. The scaling relation results 
presented here should provide a benchmark for comparison 
with cosmological simulations of the formation and evolu- 
tion of massive clusters including baryonic physics. 

Monte Carlo samples encoding the results of this pa- 
per and Paper I will be made available for download on the 
web.^'' 
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APPENDIX A: THE EFFECTS OF SELECTION 
BIASES ON THE SCALING RELATIONS 

Here we offer a pedagogical look at the effects of Malmquist 
and Eddington biases in the context of fitting the scal- 
ing relations. Consider the upper-left panel of Figure Al, 
which shows a fictitious set of cluster masses and luminosi- 
ties (crosses) and the associated luminosity-mass relation 
(red line). The simulated clusters are distributed uniformly 
in log-mass, and scattered in luminosity about the nominal 
relation following a log-normal distribution. The simplest 
way to visualize the effect of Malmquist bias on this distri- 
bution of clusters is to consider a threshold luminosity below 
which no clusters can be detected; this results in the "ob- 
served" data set shown above the dashed, blue line in the 
upper-right panel of the figure (assuming complete follow- 
up and ignoring any additional scatter due to measurement 
errors) . The distribution of observed clusters in the log-mass 
log-luminosity plane clearly does not follow the underlying 
relation indicated by the red line, and any simple method 
of fitting the data that does not account for the selection 
procedure (i.e. the luminosity threshold) will be biased by 
the low-mass clusters, which lie far from the mean relation. 

This picture is further complicated by the presence of 
Eddington bias, which results from the fact that the num- 
ber of clusters in the Universe decreases rapidly as a func- 
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Figure Al. Fictitious cluster luminosity-mass relations (red lines) and simulated data (crosses) intended to illustrate the effect of 
Malmquist and Eddington biases on the scaling relation data. In the top panels, clusters are distributed uniformly in log-mass, whereas 
in the bottom panels the distribution of log-masses is exponential. The left-hand panels reflect the true distribution of all clusters in mass 
and luminosity, while the right-hand panels show only the simulated clusters with luminosities greater than a threshold value, indicated 
by the dashed, blue lines. The figure illustrates that both the sample selection function and the underlying mass function must be taken 
into account when fitting the scaling relations. 



tion of mass. Although usually thought of in terms of the 
total number of cluster detections, the steep slope of the 
mass function has the additional effect of exaggerating the 
biased nature of the scaling relation data. This can be seen 
in the lower panels of Figure Al, which show the true and 
observed distribution of cluster masses and luminosities rel- 
ative to the true relation and the luminosity threshold when 
the distribution of cluster log-masses is exponential rather 
than uniform. Straightforwardly, the presence of many more 
low-mass than high-mass clusters results in a larger fraction 
of heavily biased (i.e. more luminous than average) data. 

In practice, the effects of selection bias on our 
luminosity-mass data are less visually obvious (cf. Figure 3) 
because the data are drawn from three flux-limited sam- 
ples, each with a different flux limit. In effect, the applica- 
ble threshold £ is a function of both redshift and the ap- 
plicable flux limit, and therefore varies from cluster to clus- 
ter. Nevertheless, the data are certainly affected by selection 
bias, as reflected in the apparent offset between the data and 
the best-fitting relation from our analysis in Figure 3. We 
stress that re-measurement of cluster luminosities through 
follow-up observations does not eliminate the effects of the 



Malmquist and Eddington biases described above because 
most of the scatter in survey-measured luminosity at fixed 
mass is due to intrinsic scatter in the scaling relation, not 
measurement errors. 

Fitting methods that do not explicitly account for the 
distribution of covariates (i.e. that assume they are uni- 
formly distributed) will naturally assign disproportionate 
importance to the lowest-mass systems, whose number is 
very sensitive to the details of the mass function in addition 
to the selection function. To properly compensate for this 
effect, it is thus crucial that an analysis incorporate infor- 
mation about both the mass function and sample selection 
function, as described in Paper I. In practice, it is necessary 
to marginalize over cosmological parameters which affect the 
mass function, particularly ag. 



APPENDIX B: PREDICTING LUMINOSITY 
MEASUREMENTS 

The distribution P{l'\z,rh, I) describes the likelihood of a 
cluster having luminosity I' measured in follow-up observa- 
tions, given that it was detected in the survey at redshift z 
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and given that the mass measured in follow-up observations 
was m (see definitions in Equation 5). In practice, the spec- 
troscopic redshift measurement error is negligible, and we 
therefore set P{z\z) = S(z — z) (where 5 is the Dirac delta 
function), i.e. z = z. With this simplification, P{l'\z,rh, I) 
can be written 

X P(mJ'\m,£,t) P{I\z,£,t), (BI) 

where P{m) is proportional to the mass function at redshift 
z; P{£,t\z,m) is the stochastic, bivariate scaling relation; 
P{rh, l'\z, m, £, t) is the likelihood of measuring mass m and 
luminosity £' given true values m, £ and t; and P{I\z,£, t) is 
the probability of detecting a cluster with redshift z, lumi- 
nosity £ and temperature t in the survey. This distribution is 
closely related to the quantity Pdot defined in Section 4.1.2 
of Paper I, the primary difference being that the information 
provided by the survey flux measurement is not incorporated 
into in the expression above. 

This paper has been typeset from a TJrX/ I^Tp^jK file prepared 
by the author. 
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